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large  object  problems  and  combination  sized  object  problems. 

Volume  I  of  this  report  Is  the  User  Manual.  The  code  execution  require¬ 
ments,  input  language  and  output  are  discussed. x„ 

Volume  II  is  the  Engineering  Manual.  The  theory  and  engineering  approxi¬ 
mations  implemented  in  the  code  are  discussed.  Modeling  criterion  are 
given. 

Volume  III  is  the  Computer  Code  Documentation  Manual.  This  manual  contain! 
extensive  software  information  of  the  code. 
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INTRODUCTION 


This  manual  contains  Instructions  for  using  the  General  Electro¬ 
magnetic  Model  for  the  Analysis  of  Complex  Systems  (GEMACS)  computer 
program.  The  program  Is  a  highly  user-oriented  general  purpose  code 
designed  for  gradual  Incorporation  of  a  variety  of  techniques  for  elec¬ 
tromagnetic  analysis  of  complex  systems.  The  user  Is  assumed  to  be  an 
experienced  electromagnetics  analyst  with  a  fair  understanding  of  applied 
linear  algebra.  The  current  version  (version  3)  of  the  code  supports  all 
of  the  functions  necessary  for  using  one  thln-wlre  and  one  surface  patch 
(Method  of  Moments  (MOM))  formalism  with  or  without  Geometrical  Theory  of 
Diffraction  (GTD)  Interactions.  GEMACS  is  implemented  in  four  sequen¬ 
tially  executable  FORTRAN  programs  (called  "modules").  The  use  of  the 
modules  In  running  GEMACS  is  described  In  section  D.  The  GEMACS  code 
uses  a  high  level  language  and  provides  flexibility  of  control  over  the 
computational  sequence  by  the  user.  Error  messages,  debug  and  trace 
options,  and  other  features  are  Included  to  aid  the  user  in  Identifying 
sources  of  fatal  errors. 

1.  Method  of  Moments  Formulations 

One  MOM  formalism  used  In  the  present  code  Includes  the  thln- 
wlre  Pockllngton  Integral  equation,  pulse  plus  sine  plus  cosine  expansion 
functions,  point  matching,  and  a  charge  redistribution  scheme  at  multiple 
wire  junctions.  This  is  the  same  formalism  as  used  In  the  AMP  (Antenna 
Modeling  Program)  code*.  The  GEMACS  code  includes  most  of  the  engineering 
features  of  the  AMP  code  such  as  loading  and  ground  plane  effects. 
However,  the  range  of  applicability  of  the  moments  technique  Is  extended 
to  objects  of  larger  electrical  size  in  the  GEMACS  code  by  using  a  solu¬ 
tion  method  for  linear  simultaneous  equations  called  BMI  (Banded  Matrix 
Iteration).  The  user  must  have  a  limited  understanding  of  the  solution 


*The  Antenna  Modeling  Program  is  a  general-purpose  thin  wire  code 
developed  and  documented  by  MB  Associates  in  Antenna  Modeling  Program- 
Engineering  Manual,  July  1972,  A0-A025890. 
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method  to  ensure  convergence  and  reasonable  efficiency.  The  method  is 
documented  in  Volume  II  (GEMACS  Engineering  Manual)  of  this  report  and  In 
references  noted  therein.  The  th in-wire  MOM  can  be  used  to  solve  general 
physical  problems  Involving  actual  wires,  wire  grid  models  of  conducting 
surfaces,  or  a  combination  of  these.  Wire  grid  modeling  is  not  yet  a 
highly  defined  process.  Modeling  guidelines  developed  in  recent  studies 
are  discussed  in  Volume  II.  The  user  must  reduce  the  physical  problem  to 
a  thln-wlre  model.  The  GEMACS  code  includes  a  highly  flexible  geometry 
processor  to  aid  in  this  task.  The  user  specifies  the  frequency,  addi¬ 
tional  features  such  as  loading  or  the  presence  of  a  ground  plane,  and 
the  excitation.  Excitation  options  currently  include  plane  or  spherical 
waves,  voltage  sources  for  antennas,  or  arbitrary  excitations  on  speci¬ 
fied  Individual  wire  segments.  Load  options  currently  include  fixed  (as 
a  function  of  frequency)  lumped  loads,  series  or  parallel  RLC  networks, 
and  finite  segment  conductivity. 

A  second  MOM  formalism  is  the  use  of  the  Magnetic  Field 
Integral  Equation  (MFIE).  Here  the  surface  current  is  expanded  in  a  set 
of  pulse  expansion  functions,  except  in  the  region  of  a  wire  connection. 
Two  orthogonal  current  directions  are  assumed  for  each  surface  patch. 
Point  matching  Is  used  at  the  patch  centers.  In  the  region  of  a  wire 
connection  to  a  patch  four  subpatches  are  generated,  and  the  continuity 
of  current  equation  at  the  center  of  the  patch  takes  into  account  the 
singular  component  due  to  the  current  flowing  from  the  wire  Into  tne 
surface.  The  MFIE  therefore  provides  a  viable  alternative  to  the  wire 
grid  modeling  approach  for  a  surface.  (See  section  C.3.g  of  the  GEMACS 
Engineering  Manual.) 

GEMACS  can  also  use  the  physical  symmetry  of  a  MOM  structure  to 
decrease  matrix  fill  time  and  matrix  equation  solution  time.  The 
symmetry  may  be  either  planar  or  rotational.  Since  only  the  physical 
symmetry  Is  used,  no  restriction  is  placed  on  the  excitation  or  loading 
of  the  structure  regarding  symmetry.  (See  section  F.2  of  the  GEMACS 
Engineering  Manual.) 
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2.  Geometrical  Theory  of  Diffraction  Formulation 

The  GTD  formalism  used  In  the  present  code  makes  It  possible  to 
model  electrically  large  objects  In  bulk  Instead  of  requiring  the 
grlddlng  of  the  object  Into  a  large  number  of  wire  segments  or  patches. 
Fields  scattered  by  the  object  are  determined  by  optic  principles:  ray 
tracing  and  reflection  coefficients.  The  GTD  extension*  to  basic  geomet¬ 
rical  optics  principles  adds  diffracted  ray  paths,  including  creeping 
waves,  and  diffraction  coefficients. 

The  present  GTD  capability  In  GEMACS  contains  the  same  computa¬ 
tional  procedures  found  in  the  BSC  (Basic  Scattering  Code),**  although 
much  has  been  added  to  BSC  to  make  It  compatible  with  the  GEMACS 
Impedance  matrix  formalism.  While  GTD  is  a  very  powerful  technique,  the 
user  must  have  some  physical  Insight  Into  the  problem  he  Is  solving  In 
order  to  interpret  correctly  the  results  obtained.  This  is  particularly 
true  when  one  sees  field  discontinuities  at  object  shadow  boundaries. 

To  use  GTD  techniques,  the  user  must  reduce  his  object  geometry 
to  a  set  of  canonical  objects:  planar  plates,  a  cylinder,  and  the 
cylinder's  two  endcaps.  These  shapes  are  described  to  GEMACS  through  the 
code's  geometry  processor.  Next,  the  field  sources  which  generate  the 
field  Incident  on  the  object,  as  well  as  the  frequency  of  the  source,  are 
specified.  Rather  than  calculate  geometry  element  currents  to  obtain 
field  pattern  data  (the  MOM  process),  the  scattered  fields  are  obtained 
directly  from  the  sources  and  geometry  by  tracing  all  geometrical  optics 
paths  from  the  sources  to  the  field  points  and  reflecting  and  diffracting 
the  waves  which  follow  these  paths  from  the  surfaces,  edges,  and  corners 
of  the  geometry  elements.  With  the  GTD  method  one  may  find  the  geometry 
surface  currents  Indirectly  from  the  total  fields  at  the  surface,  but  not 
antenna  Impedances. 


*J.  B.  Keller,  "Geometrical  Theory  of  Diffraction",  j.  Opt.  soc. 
Am.,  52:116,  1962. 

**R.  J.  Marhefka  and  W.  0.  Burnside,  "Numerical  Electromagnetic 
Code  -  Basic  Scattering  Code  Part  I:  User's  Manual,"  The  Ohio  State 
University,  ’ptembe”  979. 


3 


GTO  techniques  are  used  when  detailed  structure  elements  are 
not  Important  and  when  scattered  fields  are  the  only  quantities  desired. 
While  theoretically  It  Is  possible  to  obtain  any  level  of  detail  with 
enough  GTO  geometry  elements,  there  Is  a  practical  limit  to  the  detail 
that  can  be  obtained  because  of  the  large  number  of  second  and  third 
order  Interactions  possible  among  the  GTO  elements.  In  practice, 
therefore,  the  user  Is  forced  to  limit  either  the  number  of  GTD  elements 
or  the  number  of  physics  interactions  he  wishes  to  consider.  Guidance  on 
these  and  other  modeling  considerations  is  provided  in  Volume  II. 

A  limitation  of  the  GTO  formalism  as  implemented  in  GEMACS  is 
that  no  advantage  can  be  taken  of  symmetry.  The  full  GTD  geometry  must 
be  specified  to  GEMACS.  Furthermore,  all  GTD  physics  assumes  perfectly 
conducting  geometry  elements.  Hence  no  dielectrics  or  surface  impedance 
loadings  are  permitted.  Perfectly  conducting  ground,  however,  can  be 
accommodated  by  using  a  planar  plate. 

3.  MOM/GTD  Hybrid  Formulation 

The  presence  of  both  thin  wire  and  large  conducting  elements  on 
a  geometry  to  be  modeled  presents  a  dilemma  to  the  user  if  only  MOM  or 
GTD  techniques  are  available.  For  this  type  of  problem,  the  two  solution 
methodologies  have  been  combined  into  a  hybrid  formulation  after  the 
method  of  Thiele.*  Thin  wires  are  modeled  with  MOM  wire  segments  and 
large  conducting  objects  are  modeled  with  GTD  geometry  elements.  The 
Impedance  Interactions  normally  computed  In  the  MOM  approach  are  aug¬ 
mented  with  fields  generated  from  wire  segments,  scattered  from  the  GTD 
objects,  and  received  by  the  wire  segments. 

The  formal  mechanism  of  the  hybrid  formulation  is  transparent 
to  the  user;  the  only  extra  step  is  to  divide  his  geometry  into  MOM  and 
GTD  elements.  The  same  commands  used  to  generate  the  interaction  and 
excitation  matrices  for  the  MOM  problem  are  used  for  the  MOM/GTD 
solution. 


*G.  A.  Thiele  and  T.  H.  Newhouse,  "A  Hybrid  Technique  for  Combining 
Moment  Methods  with  the  Geometrical  Theory  of  Diffraction,"  ieee  Trans¬ 
actions  on  Antennas  and  Propagation,  January  1975. 


The  primary  computational  difference  in  the  MQM/GTD  technique 
Is  the  small  size  of  the  matrices.  For  example,  the  interaction  matrix 
size  Is  determined  from  the  number  of  MOM  elements  regardless  of  the 
number,  size,  or  complexity  of  the  GTD  geometry.  Hence,  a  blade  antenna 
mounted  on  a  B-52  aircraft  would  require  the  same  size  interaction  matrix 
as  the  same  antenna  in  free  space.  Each  matrix  element,  however, 
consists  of  an  interaction  made  up  of  a  direct  (MOM)  interaction  and 
(perhaps  hundreds)  of  indirect  (GTD)  interactions.  While  each  matrix 
element  takes  longer  to  generate,  the  small  matrix  si2e  makes  computation 
of  the  solution  more  efficient  than  attempting  to  grid  the  entire 
geometry  with  MOM  objects  and  solving  a  huge  set  of  simultaneous 
equations. 

The  restrictions  on  the  hybrid  formulation  contain  all  the 
restrictions  of  both  MOM  and  GTD  formulations.  The  user  is  referred  to 
Volume  II  for  a  discussion  of  these  and  the  few  additional  restrictions 
applicable  to  the  hybrid  method. 

4.  Numerical  Solution  Methods 

When  executing  a  MOM  or  MOM/GTD  problem,  the  code  generates  a 
set  of  linear  simultaneous  equations  from  the  information  provided.  The 
user  controls  the  process  by  which  the  equations  are  solved.  If  the 
total  number  of  wire  segments  in  the  model  is  sufficiently  small, 
standard  solution  methods  are  applicable.  Solution  by  full  matrix 
triangular  decomposition  is  one  of  the  least  expensive  general  methods 
and  is  supported  by  GEMACS.  For  large  problems,  even  this  method  is  too 
expensive,  and  the  BMI  solution  method  should  be  specified  by  the  user. 
This  method  is  considerably  less  expensive  provided  the  user  carefully 
chooses  the  segment  numbering  and  matrix  bandwidth  according  to  the 
guidelines  discussed  in  Volume  II. 

The  user  may  specify  other  quantities  to  be  computed  from  the 
wire  currents,  such  as  impedances,  coupling  parameters,  near  fields  or 
far  fields.  These  are  computed  from  currents  regardless  of  the  solution 
process  specified.  Regardless  of  the  solution  technique  exercised,  it  is 
emphasized  that  the  user  must  be  familiar  with  general  results  from  the 


literature  to  ensure  that  the  computed  solution  using  the  model  for  the 
system  is  of  sufficient  accuracy  for  the  intended  purpose.  For  example, 
the  far  fields  can  be  computed  from  approximate  currents  obtained  by 
specifying  a  weak  convergence  criterion  when  using  the  BMI  solution 
method.  This  will  allow  the  reduction  of  the  required  computer  resources 
when  large  systems  are  being  analyzed.  Near-field  computations  will 
require  more  accurate  wire  currents,  obtainable  only  with  stricter,  more 
costly  convergence  criteria. 

The  present  code  generates  an  interaction  matrix  from  the  MFIE 
and  EFIE  (Electric  Field  Integral  Equation)  discussed  in  the  GEMACS 
Engineering  Manual.  The  wire  current  is  represented  by  a  sine,  cosine, 
and  pulse  expansion  function  with  redistribution  at  junctions  based  on 
the  fractional  length  of  each  segment  with  respect  to  the  total  length  of 
all  segments  connected  at  the  junction.  The  surface  current  is  repre¬ 
sented  by  a  pulse  function.  The  interaction  matrix  may  be  modified  by 
loading  the  individual  wire  segments  of  the  model  using  resistance, 
capacitance,  and  inductance  in  parallel  or  series  configurations. 

Associated  with  the  geometric  structure  and  interaction  matrix 
is  an  excitation  matrix  which  contains  the  total  tangential  electric 
field  present  at  the  midpoint  of  each  segment  or  patch.  The  electric 
field  may  be  caused  by  as  many  combinations  of  three  types  of  sources  as 
desired.  These  types  are  plane  and  spherical  wave  sources  for  scattering 
problems  and  voltage  sources  for  antenna  problems.  In  addition,  the  user 
may  assign  an  arbitrary  value  to  the  excitation  of  any  wire  segment  to 
force  the  desired  boundary  condition. 

With  the  interaction  matrix  denoted  by  [ZJ  and  the  excitation 
matrix  denoted  by  (El,  the  primary  function  of  the  code  is  to  generate 
and  solve  the  system  of  equations  for  the  electric  current  [I]. 

[Z]  [I]  =  [E] 


This  may  be  done  using  direct  full  matrix  decomposition  in  which  [ Z 1  is 
decomposed  into  lower  and  upper  triangular  matrices  [ ZL ]  and  [ ZU ] .  For¬ 
ward  elimination  and  back  substitution  are  then  performed  as  indicated 
below. 


o 

.  \  i 

.*  v.\ 


.  *  » 


***  •/*/  •//.*//  .•*/.%*.  V  >  / 
■  V* .  *Y  V*V-  \  • 


■*(* 


[Z  ]  [I  ]  =  [E  ] 

[ZL]  [ZU]  [I  ]  =  [E  ] 

[ZLJ  [I*]  =  [E  1 

[ZU]  [I  ]  =  [I ' ] 


where,  since  [ZL]  and  [ZU]  are  triangular,  the  actual  Inverse  Is  not 
required. 

For  very  large  problems,  the  direct  solution  method  may  be 
prohibitive  due  to  the  large  amount  of  time  required  and  the  possible 
roundoff  errors.  In  this  case,  the  Banded  Matrix  Iteration  (BMI)  tech¬ 
nique  Is  available.  In  this  method,  [Z]  is  partitioned  (not  decomposed) 
Into  lower  and  upper  triangular  matrices  ([LZ]  and  [UZ] )  and  a  central  or 
banded  matrix  [BZ] .  This  Is  Illustrated  graphically  below. 


[Z] 


[Z]  «  [LZ]  +  [BZ]  +  [UZ] 


where  LZ',  BZ',  UZ'  are  the  nonzero  elements  of  LZ,  BZ,  and  UZ.  In  this 
notation,  the  problem  to  be  solved  becomes: 


or: 


[Z] 

[I] 

-  [LZ] 

[I]  + 

[BZ]  [I]  +  [UZ]  [I]  =  [E] 

[BZ] 

[I] 

=  [El 

-  ([LZ] 

+  [UZ])  [I] 

The  banded  portion  of  Interaction  matrix  [BZ]  should  contain  the 
dominant  Interactions  while  those  contained  In  [LZ]  and  [UZ]  may  be 
viewed  as  perturbations  In  a  properly  posed  problem.  In  this  case,  the 
BMI  technique  Involves  appending  subscripts  to  the  [I]  vector  and  only 
solving  the  banded  portion  of  the  Interaction  matrix. 
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[BZ]  [II 


[RHS1 


where  the  right-hand  side  at  the  n  Iteration  [RHS] n  is  given  by: 

[ RHS ] n  =  [t]  -  ( [LZ]  +  [UZJ)  [I]nl 

The  starting  value  for  [I]  or  E 1 1 0  Is  zero  unless  preset  by  the  user. 
There  is  very  little  if  any  advantage  to  presetting  [I].  Since  [BZ]  Is 
usually  much  smaller  than  [Z],  the  time  to  perform  lower/upper  decompo¬ 
sition  is  reduced  considerably  and  the  system  to  be  solved  becomes: 


[BZL]  [BZU]  [I]n  *  [RHS]n 


where  [BZL]  and  [BZU]  are  the  lower  and  upper  triangular  matrices 
obtained  by  decomposing  [BZ] .  The  solution  for  [ I ]  n  is  obtained  exactly 
as  for  the  full  matrix  decomposition.  When  using  BMI,  the  user  must 
provide  the  convergence  measure  and  value  to  be  used  to  stop  the  itera¬ 
tive  procedure.  Three  criteria  or  measures  are  available,  the  BCRE 
(Boundary  Condition  Relative  Error),  the  IRE  (Iterative  Relative  Error), 
and  the  PRE  (Predicted  Relative  Error).  The  BCRE  Is  a  measure  of  how 
well  the  solution  matches  the  boundary  condition.  Mathematically: 


BCRE  * 


I  EE]  -  [Z]  [iy 


While  this  form  has  great  engineering  appeal,  mathematically  It  Is  not 
recommended  since  the  system  of  equations  may  be  ill-conditioned  near 
resonances,  and  there  Is,  by  definition,  a  large  variance  In  the  elements 
of  [I]  which  will  result  In  a  small  BCRE. 

The  second  criterion  Is  the  IRE.  This  Is  defined  mathemat¬ 
ically  as: 


'<•1,  -  O'n-l1 


and  can  be  seen  to  be  the  relative  change  between  successive  approxima¬ 
tions  to  the  solution.  For  slowly  converging  problems,  this  criterion 
may  cause  premature  termination  of  the  Iterative  procedure. 

Finally,  the  PRE  may  be  used.  This  quantity  Is  determined  by 
using  an  exponential  fit  to  the  two  previous  values  for  the  IRE  and  has 
been  shown  to  be  a  good  approximation  to  the  ARE  (Actual  Relative  Error) 
after  four  Iterations.  The  ARE  Is  defined  as: 


ARE 


i iH,  -  mi 

mil 


where  [I]  Is  the  exact  solution.  Since  the  exact  solution  is  not  avail¬ 
able,  the  PRE  is  the  recommended  criterion.  (See  the  GEMACS  Engineering 
Manual,  section  F.1.3  for  a  discussion  of  the  PRE.) 

The  value  of  the  convergence  criterion  depends  largely  on  the 
output  desired.  If  Input  Impedance  or  near-field  parameters  are  desired, 
a  1  percent  value  Is  not  Inappropriate;  however.  If  normalized  far-fleld 
patterns  are  desired,  a  10  to  15  percent  value  may  be  sufficient. 

Once  the  solution  has  been  obtained,  the  Input  Impedance  of 
each  voltage  driven  element  (l.e.^  Antenna  Feed  Point)  Is  output  to  the 
user.  These  are  computed  as  Za  =  -p  since  a  delta  gap  model  Is  used  for 
antenna  sources.  The  currents  may  ^>e  used  as  inputs  to  the  field  compu¬ 
tation  routines  to  obtain  the  near-  and/or  far-electric  field  patterns, 
and  the  coupling  between  pairs  of  antennas.  The  user  should  note  that 
the  Impedance  value  printed  Is  equal  to  the  antenna  Input  Impedance  only 
when  a  single  voltage  driven  element  Is  specified.  Otherwise  Z.  Is 
merely  the  ratio  of  element  voltage  to  element  current  and  cannot  be 
related  to  an  Input  Impedance. 

The  above  discussion  Is  relevant  only  to  MOM  or  MOM/GTD  hybrid 
problems.  A  problem  utilizing  only  GTD  geometries  does  not  require  a 
numerical  solution,  since  no  simultaneous  equations  are  generated. 
Rather,  the  GTO  solution  (radiated  and  scattered  fields)  is  obtained 
directly  from  the  ray  tracing,  reflection,  diffraction,  and  shadowing 
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computations  that  take  place  where  GTO  fields  are  requested.  In  this 
case  surface  currents  and  Input  Impedances  are  not  calculated  by  the 
code. 

5.  GEHACS  Capabilities  Summary 

The  present  version  of  GEHACS  supports  three  electromagnetic 
solution  methods  and  two  numerical  techniques.  These  are: 

Method  of  Moments  (wire  grlddlng,  patches) 

Geometrical  Theory  of  Diffraction  (Ray  tracing,  optics, 
diffraction) 

A  MOM/GTD  Hybrid 

Matrix  Triangular  Decomposition 

Banded  Matrix  Iteration 

There  are  Inherent  limitations  to  the  solution  techniques  available.  The 
user  who  Is  not  familiar  with  these  techniques  should  consult  the  Engi¬ 
neering  Manual  and  Its  references  to  avoid  wasting  valuable  time  and 
computer  resources  working  an  111-posed  problem. 


mu# 


B.  COMPUTATIONAL  APPROACH 


The  basic  approach  In  the  design  of  the  GEMACS  code  Is  to  permit  the 
user  to  generate  or  define  data  sets  and  then  to  perform  operations  on 
the  data  contained  In  the  data  sets.  The  data  sets  are  Identified  using 
symbolic  names  of  six  ANSI  FORTRAN  characters  or  less.  All  names  must 
start  with  an  alphabetic  character.  The  result  of  this  approach  is  a 
code  which  the  user  completely  controls  down  to  the  functional  or  opera¬ 
tional  level.  Associated  with  each  symbol  Is  a  set  of  symbol  character¬ 
istics  referred  to  as  attributes.  These  attributes  are  generated  as  the 
data  associated  with  the  symbol  are  generated  or  modified  by  an  opera¬ 
tion.  The  attributes  are  checked  each  time  the  data  are  used  in  an 
operation.  This  Insures  the  integrity  of  the  resulting  data  and  the 
sanity  of  the  operation.  For  example.  If  the  user  attempts  to  generate 
an  Interaction  matrix  using  data  which  do  not  represent  a  geometrical 
structure,  an  error  will  occur  since  the  operation  is  not  defined. 
However,  If  an  Interaction  matrix  is  generated  for  a  properly  defined 
geometry  data  set.  It  will  have  the  attributes  of  a  complex  interaction 
matrix  and  will  be  Identified  as  having  been  generated  from  the  geometry 
data  set  specified.  In  this  way,  a  symbol's  lineage  and  the  type  of  data 
(both  physical  and  computational)  are  known.  Likewise,  the  solution 
vector  for  a  geometry  data  set  is  linked  to  the  interaction  matrix  and 
thus  to  the  geometry  which  was  the  parent  data  for  all  subsequent  opera¬ 
tions.  Therefore,  when  field  output  Is  desired,  the  code  will  retrieve 
the  correct  electrical  current  data  for  use  with  the  geometry  data 
specified. 

The  primary  function  of  the  GEMACS  overhead  routines  is  to  store  and 
retrieve  the  data  sets  In  order  that  user  specified  operations  may  be 
performed.  The  result  of  these  operations  will  be  a  solution  or  analysis 
of  the  system  described  to  the  code. 

As  more  operations  are  added  to  the  GEMACS  code,  it  is  necessary  to 
define  the  attributes  of  the  data  needed  to  assure  correct  program  opera¬ 
tion,  modify  the  Input  Module  to  recognize  the  command  and  add  the 


software  to  perform  the  operation.  The  data  handling  will  be  taken  care 
of  by  the  executive  level  programs  in  the  GEMACS  code. 

As  an  example,  a  normal  MOW  scattering  solution  at  a  single 
frequency  for  which  the  far-field  radiation  pattern  is  desired  could  be 
done  with  the  following  input  stream.  The  directives  are  discussed  in 
the  following  section. 

COMMAND 

1  GMOATA  =  GE0M1 

2  EINC  *  ESRC  (GE0M1),  FRQ  =  180.0,  SW  =  l.,0., THETA  =  45. 

3  ZGEN  GMDATA  =  GE0M1,  ZMATRX  *  ZMAT1 

4  BNDZl  =  BAND  (ZMAT1) ,  BNDW  *  50 

5  BNDZl  =  LUO  (BNDZl) 

6  BNDZl  *  CUR1  «  EINC  -  ZMAT1  *  CUR1,  MAXITR  =  10,  CONVRG  *  PRE, 
VALUE  =  20 

7  EFIELO  (CUR1) ,  LOGPLR,  PI  =0.  P2  -  180.  DP  *  10.  T1  *  90. 

8  END 
{GEOMETRY  DATA} 

END 

Card  1  directs  the  geometry  processor  to  generate  the  data  to  be  asso¬ 
ciated  with  symbol  GE0M1.  Card  2  directs  that  data  associated  with 
symbol  EINC  be  generated  by  a  spherical  wave  at  180  MHz  Incident  on  the 
geometry  specified  by  GE0M1  at  a  spherical  angle  theta  of  45°.  Card  3 
directs  that  data  associated  with  Interaction  matrix  ZMAT1  be  generated 
using  the  sine  +  cosine  +  pulse  expansion  function  on  the  geometry  asso¬ 
ciated  with  GE0M1.  Card  4  causes  extraction  of  the  elements  from  ZMAT1 
which  are  located  within  50  elements  of  the  main  diagonal  and  associates 
the  data  with  symbol  BNDZl.  Card  5  results  In  BNDZl  being  decomposed 
Into  upper  and  lower  triangular  matrices  and  the  result  stored  In  BNDZl. 


Card  6  Invokes  the  BMI  solution  technique  to  obtain  the  data  for  symbol 
CUR1.  The  procedure  Is  limited  to  10  Iterations  and  the  PRE  convergence 
criterion  will  be  used.  When  the  PRE  <  20  percent,  the  procedure  will 
stop  and  return  to  execute  the  next  directive.  Card  7  directs  the 
computation  and  output  of  the  far  field  of  structure  GE0M1.  The  field 
will  be  computed  at  10°  Intervals  from  t  =  0  to  $  *  180  for  0  *  90°.  In 
addition  to  a  tabular  print,  the  logjQ  of  the  field  will  be  plotted  on  a 
polar  graph.  Card  8  Indicates  the  end  of  the  directives  and  Is  followed 
by  the  geometry  data  Input  which  Is  also  terminated  by  an  END  card.  The 
Input  Is  completely  free  field  and  there  are  default  values  for  most 
parameters. 

All  of  the  directives  are  processed  before  any  execution  begins. 
This  precludes  wasting  considerable  computer  resources  In  the  event  of  an 
error  In  the  command  directives. 


C.  GEMACS  COMMANDS  AND  GEOMETRY  LANGUAGE 


The  GENACS  Inputs  are  In  two  categories.  The  command  language 
directs  the  program  execution  while  the  geometry  language  describes  the 
geometrical  properties  of  the  structure  being  analyzed. 

The  GEMACS  command  language  Is  a  free-fleld,  keyword  oriented  input 
stream.  The  order  of  the  data  Is  generally  not  Important,  and  the  items 
on  each  card  are  delimited  by  a  blank  or  a  comma.  An  Item  Is  considered 
to  be  all  of  the  Input  associated  with  a  particular  parameter  such  as 
THETA  *  90.  Note  that  an  Item  may  consist  of  several  entries,  and  each 
entry  Is  referred  to  as  a  field.  Blanks  may  be  imbedded  between  fields 
of  an  item  but  not  within  a  field.  Thus,  THETA  =  90.  Is  acceptable  while 
THETA  *  9  0.  Is  Incorrect  and  will  be  Interpreted  as  two  Items  (i.e., 
THETA  *  9  and  0).  This  extraneous  Item  would  be  detected  by  the  code  and 
execution  would  be  Inhibited. 

In  order  to  prevent  wasting  computer  resources,  all  of  the  GEMACS 
commands  are  read  by  the  INPUT  module  and  checked  for  errors  prior  to 
execution  of  the  code.  All  errors  in  the  user  commands  are  Identified; 
that  is,  one  error  does  not  preclude  location  of  any  other  error  during 
the  same  execution.  This  prevents  the  need  for  the  user  to  make  several 
submittals  to  debug  the  Input. 

The  GEMACS  geometry  language  is  also  a  free-fleld  language. 
However,  the  items  must  appear  In  the  order  specified  or  an  error  will 
occur  which  may  not  be  detected.  The  reason  for  not  using  keyword- 
specified  Items  on  the  geometry  Inputs  Is  to  decrease  the  effort  required 
by  the  user  since  the  geometry  Inputs  are  usually  much  larger  than  the 
command  Inputs. 

For  both  input  types  there  are  several  standard  conventions.  These 
relate  to  comment  cards,  comments  on  cards,  and  continuation  cards. 


Comment  cards  are  those  cards  which  have  a  $  as  the  first  nonblank 
character.  Likewise,  comments  may  be  appended  to  a  command  or  geometry 
Input  by  preceding  the  comment  with  a  $•  If  the  last  character  encoun¬ 
tered  before  a  $  or  the  end  of  a  card  is  a  comma  or  arithmetic  operator 
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(+,-,*,/),  the  next  card  must  be  a  continuation  card.  If  a  card  has  a 
continuation  character  In  column  1,  It  Is  treated  as  a  continuation  of 
the  previous  card.  A]2  continuation  cards  must  have  a  continuation 
character  In  column  1.  The  continuation  character  may  be  Installation 
dependent.  It  Is  an  asterisk  (*)  In  version  3.  Other  possible  choices 
are  the  other  arithmetic  operators.  The  actual  character  is  defined  as 
variable  NCONCH  and  Is  set  In  subroutine  BLKDAT. 

1.  GEMACS  Command  Language 

In  describing  the  GEMACS  command  language.  Items  enclosed  In 
brackets  [  ]  have  default  values  and  need  not  be  specified  If  the  default 
value  Is  acceptable  to  the  user.  Items  enclosed  In  braces  {  }  indicate  a 
multiple  choice.  The  only  restrictions  on  symbolic  names  provided  by  the 
user  are  that  they  be  six  characters  or  less,  the  first  character  be  a 
member  of  the  alphabet  (A-Z),  and  only  characters  contained  In  the  alpha¬ 
bet  or  the  digits  (0-9)  be  present.  That  Is,  the  characters  =,  +,  -,  *, 
/,  $,  and  comma  are  not  allowed.  In  addition,  the  following  reserved 
keywords  may  not  be  used  for  symbolic  names. 

ONE  LETTER  KEYWORDS 


A 

C 

0 

N 

0 

R 

V 

X 

z 

TWO  LETTER  KEYWORDS 

CO 

CR 

CS 

CW 

CY 

Cl 

C2 

DC 

DM  DP 

DR 

OT 

OW 

OX 

OY 

DZ 

EC 

ED 

El  ER 

ES 

EU 

IS 

LU 

MM 

NP 

NR 

ON 

PC  PD 

PL 

PR 

PI 

P2 

RC 

RO 

RR 

R1 

R2  SC 

SW 

TI 

T2 

VS 

XI 

X2 

Y1 

Y2 

Z1  Z1 

THREE  LETTER 

KEYWORDS 

ABS 

COP 

ECC  ENO 

FRQ 

GTD 

ILP 

INV 

LUO 

OFF 

POR  PHI 

PRE 

RDP 

SEQ 

SET 

FOUR 

LETTER  1 

KEYWORDS 

AXIS 

BAND 

BCRE 

BNDW 

CONO 

EPSR 

ESRC 

LOOP 

PLOT 

PRLC 

READ 

SCOP 

SEGS 

SIZE 

SRDP 

SRLC 

TAGS 

TIME 

TYPE 

VSRC 

ZGEN 

ZIMP 
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FIVE  LETTER  KEYWORDS 


CONJG 

CPINC 

CPNUM 

OEBUG 

INPUT 

LABEL  PARTN 

PIVOT 

PRINT 

PULSE 

PURGE 

SOLVE 

STATS 

THETA 

TRACE 

VALUE 

WRITE 

SIX  LETTER  KEYWORDS 


8ACSUB 

CHKPNT 

COLPSE 

CONVRG 

EFIELD 

EXPAND 

FILE  ID 

GMOATA 

LINLIN 

LINLOG 

LINPLR 

LOGLIN 

LOGLOG 

LOGPLR 

MAXITR 

MODULE 

NUMFIL 

OUTPUT 

PCESIN 

REDUCE 

REFLCT 

REPLAC 

RSTART 

SETINT 

SINCOS 

SYMDEF 

TRANSP 

WIPOUT 

ZCODES 

ZLOADS 

ZMATRX 

addition 

to  these 

names , 

three 

other  cases  are 

avoided.  If  an  out-of-core  matrix  Is  to  be  decomposed  and  is  identified 
symbolically,  as  XXXXXX  where  the  Xs  are  legitimate  characters,  then  the 
user  must  not  define  another  data  set  with  either  XXXLWR,  XXXUPR,  or 
XXXPVT.  These  names  will  be  Internally  generated  to  contain  the  lower 
triangular  decomposed  matrix,  the  upper  triangular  decomposed  matrix,  and 
the  pivot  vector  as  specified.  These  names  may  be  referenced  in  output 
statements.  However,  In  general,  the  user  should  simply  ignore  them.  It 
Is  also  correct  to  assume  that  a  matrix  which  resides  out  of  core  is  no\ 
destroyed  when  It  Is  decomposed.  However,  unless  a  matrix  is  large,  it 
will  be  difficult  for  the  user  to  know  a  priori  where  it  will  be  stored. 

The  second  case  Is  when  an  Interaction  matrix  is  generated  with 
GTD  Interactions  specified.  A  shadowing  matrix  Is  generated  by  the  code, 
whose  symbolic  name  Is  XXXSHO,  where  the  Xs  are  the  leftmost  characters 
of  the  Interaction  data  set  name. 

The  third  case  involves  generation  of  GTD  scattered  fields. 
Here,  the  code  generates  a  Green's  function  scattering  matrix  XXXGFM, 
where  the  Xs  are  the  leftmost  characters  of  the  field  matrix  data  set 
name. 

A  list  of  commands  is  given  in  table  1  In  which  the  following 
codes  are  used: 

S  -  Previously  undefined  symbol 
DS  -  Previously  defined  symbol 
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Either  S  or  OS 
Numeric  value 

Either  a  DS  or  numeric  value 
Keyword 

Alphameric  word  of  at  most  six  characters  starting  with  a 
character  A-Z  and  containing  only  characters  A-Z  and  0-9. 


COMMANDS 


COMMANDS,  FORMATS,  AND  MNEMONICS  (Continued) 


COMMANDS,  FORMATS,  AND  MNEMONICS  (Concluded) 


V  , 
U3I 


BACSUB 


Forward  EHminatlon/Back  Substitution  Command 
BACSUB  DS1*SDS  =  DS2 


This  command  causes  the  solution  vector  for  a  previously 
decomposed  matrix  0S1  to  be  found  for  the  excitation  vector  DS2  and 
stored  as  symbol  SOS.  For  example,  assuming  SRC1,  SRC2,  and  SRC3  have 
been  previously  defined,  then  the  solutions  SOLI,  S0L2,  and  S0L3  would  be 
obtained  by: 

ZM  =  LUD(ZIJMAT) 

BACSUB  ZM*S0L 1  =  SRC1 
BACSUB  ZM*S0L2  =  SRC2 
BACSUB  ZM*S0L3  =  SRC3 

This  command  permits  the  user  to  obtain  several  solutions  with  only  one 
decomposition.  Note  that  this  is  in  contrast  to  multiple  uses  of  the 
SOLVE  command,  each  of  which  would  initiate  its  own  decomposition 
process. 
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Construct  Banded  Matrix 


SOS  =  BAND  (DS),  BNDW  =  N 


This  operation  causes  the  data  associated  with  the  matrix  DS 
which  are  within  N  elements  of  a  diagonal  element  to  be  transferred  to 
the  symbol  identified  as  SDS.  Note  that  SDS  and  DS  may  not  be  the  same 
symbolic  name.  This  operation  is  typically  used  to  construct  the  banded 
matrix  for  use  in  the  BMI  solution  process.  This  operation  is  illus¬ 
trated  in  the  figure. 

Examples: 

BNDZIJ  =  BAND  (ZIJMAT),  BNDW  =  50 

This  operation  will  construct  a  banded  matrix  from  the  data 
associated  with  ZIJMAT. 


DIAG  *  BAND  (ZIJMAT),  BNDW  *  0 

This  operation  will  extract  the  diagonal  elements  from  ZIJMAT 
and  store  them  as  DIAG. 


BNDW+] 


I 

\lN 


\  SDS 
\ 

S 

\ 


Illustration  of  BAND  Operation 


BMI  Solution  Process 


IBCRE  ) 

0S1  *  S0S1  =  0S2  -  DS3  *  SDS2,  MAXITR  =  N  .CONVRG  =  IRE 

I  p«e  t 


[.VALUE  =  N  J 

This  command  executes  the  BMI  solution  process.  DS1  must  be  a 
banded  decomposed  matrix  whose  elements  were  originally  contained  In  0S3. 
The  solution  will  be  stored  In  S0S1  upon  completion.  0S2  Is  the  excita¬ 
tion  or  right-hand  side  of  the  original  system  of  simultaneous  equations 
and  DS3  Is  the  original  Interaction  matrix  of  coefficients.  DS3  will 
still  contain  the  elements  which  were  used  to  generate  the  banded  matrix, 
but  they  will  be  Ignored.  The  symbol  S0S2  may  be  predefined  and 
preloaded  or  It  may  be  S0S1,  In  which  case  It  will  be  Initialized  to 
zero.  The  convergence  parameter  to  be  used  Is  contained  In  the  CONVRG 
Item.  Their  relative  merits  and  definitions  of  the  three  convergence 
parameters  are  described  In  Volume  II,  GEMACS  Engineering  Manual. 
Default  convergence  Item  is  CONVRG  =  PRE.  The  value  In  percent  which  the 
convergence  parameter  must  reach  Is  contained  In  the  VALUE  Item.  The 
default  Item  Is  VALUE  =1.  The  MAXITR  parameter  defines  an  upper  bound 
on  the  number  of  Iterations.  There  Is  no  default  value  for  this  Item. 
It  must  appear  before  the  CONVRG  and  VALUE  Items. 

Example; 

LUDZIJ  *  CUR  =  VDRV  -  ZIJMAT  *  CUR,  MAXITR  =  5  CONVRG  =  BCRE, 


CHKPNT 


Checkpoint  Command 
CHKPNT  [LU  =  N]  [.CPINC  =  N]  [  ,NR] 

This  command  designates  the  FORTRAN  logical  unit  (LU)  N  In  the 
LU  *  N  Item  to  be  used  to  receive  the  checkpoint  data.  The  default  Item 
Is  LU  *  7,  and.  If  the  user  specifies  a  different  logical  unit,  he  must 
assure  the  availability  of  the  unit  to  the  GEMACS  code.  The  first  word 
written  on  the  checkpoint  file  will  be  the  name  of  the  module  generating 
the  checkpoint.  Checkpoints  will  be  taken  at  time  increments  of  N  CP 
minutes  specified  in  the  CPINC  =  N  item.  If  the  item  is  not  specified, 
an  Immediate  checkpoint  is  written  when  the  command  is  encountered  during 
execution.  This  type  of  command  will  not  change  the  checkpoint  increment 
specified  on  a  previous  command.  Multiple  CHKPNT  commands  may  be  used  to 
vary  the  checkpoint  increment  and  logical  unit  during  execution. 

Checkpoints  are  accomplished  by  writing  all  data  contained  in 
named  commons  and  all  data  associated  with  symbolic  names  to  the  check¬ 
point  file.  Large  problems  can  generate  a  very  large  amount  of  data  and 
it  is  advisable  to  avoid  using  magnetic  tapes  for  the  checkpoint  file 
since  multiply  reels  may  be  required.  Also,  a  large  CP  increment  is 
recommended  unless  large  data  blocks  are  PURGED  when  no  longer  needed. 

A  historical  record  of  checkpoint  information  is  kept  if  the  NR 
parameter  is  specified;  otherwise,  the  checkpoint  file  is  rewound  after 
each  checkpoint  and  overwritten  with  subsequent  checkpoints.  Due  to  the 
large  amount  of  data,  use  of  the  NR  parameter  is  not  recommended.  If  the 
NR  parameter  is  used,  checkpointing  should  be  controlled  directly  from 
the  command  language  by  omission  of  the  CPINC  item.  This  is  the  only 
mode  In  which  restarting  can  be  guaranteed  with  known  data  for  multiple 
checkpoints  on  the  same  file  in  version  3  of  GEMACS. 

The  default  checkpoint  file  is  rewound  prior  to  writing  the 
end-of-module  checkpoint.  If  intermediate  checkpoints  are  to  be  saved, 
an  alternate  logical  unit  must  be  specified  in  the  LU  =  N  item. 

25 


.v.  >•.  *■.  v-v-  V-  v.v  w. - 

.  *  *  ,  .  »  .  • .  v  Jr*  A  .  •  . » 


DEBUG 


Debug  Command 

SON  j 

OFF  (  [ ,  IIP] 

TRACE  ? 

STATS) 

This  command  is  used  to  obtain  diagnostic  information  during 
program  execution.  Specifying  the  ON  parameter  causes  all  available 
Information  associated  with  the  subsequent  tasks  to  be  printed.  When  the 
TRACE  parameter  is  specified,  the  printout  will  include  subroutine  entry 
and  exit  information  to  allow  the  user  to  follow  the  program  flow.  When 
STATS  is  specified,  the  subroutine  timing  and  entry  statistics  are  accum¬ 
ulated  and  output  upon  completion  of  the  execution.  When  OFF  is 
specified,  the  program  returns  to  the  normal  mode.  The  parameter  IIP  may 
be  specified  if  the  user  needs  to  obtain  diagnostic  information  during 
the  execution  of  the  Input  Module. 

Examples: 

1.  DEBUG  ON 

BACSUB  Z  *  I  =  V 
DEBUG  OFF 

This  command  stream  will  cause  a  detailed  printout  to  occur 
during  execution  of  the  BACSUB  command. 

2.  DEBUG  ON, ILP 

BACSUB  Z  *  I  =  V 
DEBUG  OFF 

This  will  cause  a  detailed  printout  to  occur  during  input 
processing  and  during  execution  of  the  BACSUB  command. 

3.  DEBUG  ON, ILP 
END 

This  will  cause  all  of  the  input  language  processes  and  execu¬ 
tion  tables  to  be  printed  on  termination  of  input  processing  and  before 
execution  of  the  tasks  specified  by  the  user. 


•vrv 


a 
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Arithmetic  Operation  (Scalar  Quantities  Onl 


SOS  =  0SN1 


This  command  directs  the  arithmetic  operation  specified  to  be 
performed  on  the  data  associated  with  the  symbolic  name  or  the  numeric 
value  used.  Note  that  the  resultant  symbol  may  be  the  same  as  an 
operand.  Several  global  internal  parameters  may  be  defined  by  use  of 
arithmetic  operations.  These  are: 

FRQ  (frequency  in  MHz) 

TIME  (CP  run  time  in  minutes) 

NUMFIL  (highest  FORTRAN  logical  unit  number  available  for 
use) 

COND  (ground  conductivity  (in  mhos/m)) 

EPSR  (relative  dielectric  constant  for  ground) 

In  version  3  of  GEMACS  there  is  no  hierarchy  of  operation  and 
the  user  must  use  parentheses  to  denote  order  of  operations.  Operations 
are  performed  from  right  to  left  and  from  innermost  to  outermost  paren¬ 
thesis.  Thus, 


while 


A  =  3**2  +6+3 

A  =  (3**2)  +  6  +  32  +  6 

Examples: 

FRQ  =  FRQ  +  1 . 

OMEGA  =  6.28  *  FRQ 
VOLTS  =  ZIJMAT  *  CURENT 
TIME  =  30 


SINCREMENT  FREQUENCY 
{CONVERT  TO  RADIAN 
{MULTIPLY  TWO  SCALARS 
{SET  30  MINUTE  CP  TIME  LIMIT 


. w  .v mv.v  -n  .vy.j.%  y>v,v.' 


*-vV' 

->  -  - . 
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EFIELD 

Electric  Field  Output 


SOS  =  EFIELD  (OS)  , 


[,U2=N]  [ ,DU=N]  [,V2=N]  [,DV=N]  [,W2=N] 

[ ,DW=N]  ( ,U1=N]  [ ,V1=N]  [ ,W1=N] 

This  command  will  compute  the  incident,  scattered,  or  total 
electric  field  due  to  the  data  set  identified  as  DS.  The  resultant  data 
are  associated  with  the  symbol  In  the  SDS  Item. 

The  field  quantity  computed  is  determined  by  the  data  set  type 
of  OS  (solution  or  source)  and  the  Interactions  set  by  SETINT.  Normally 
the  scattered  field  Is  computed.  However,  if  El,  ES,  or  EU  has  been 
specified  in  the  last  SETINT  command,  incident  or  total  fields  may  be 
calculated.  In  this  case,  total  fields  (incident  plus  scattered  fields) 
are  calculated  If  DS  is  a  solution  data  set;  if  DS  is  a  source  data  set 
only  incident  fields  are  computed. 

The  location  of  all  points  at  which  the  field  is  to  be  computed 
mey  be  specified  In  spherical,  cylindrical,  or  Cartesian  coordinates.  If 
spherical  parameters  are  specified  and  R  is  omitted,  the  far  electric 
field  will  be  computed.  If  R  Is  specified,  the  near  field  will  be 
computed.  The  order  in  which  the  parameters  are  specified  will  determine 
the  order  of  the  output.  Ul,  VI,  and  W1  may  be  Rl,  Tl,  PI,  XI,  Yl,  Zl, 
with  1)2,  V2,  and  W2  being  R2,  T2,  P2,  X2,  Y2,  Z2;  and  DU,  DV,  and  DW 
being  DR,  DT,  DP,  DX,  DY,  DZ.  In  this  way,  the  user  specifies  the  first 
point,  the  increment,  and  the  last  point  for  each  coordinate  axis.  If  U, 


LINLIN 

LINLOG 

LOGLIN 

LOGLOG 

LINPLR 

LOGPLR 
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EFIELD  (Continued) 


W,  V  represent  three  coordinate  specifications,  then  specifying  U1  before 
W1  will  cause  the  variation  specified  for  W  to  be  completed  for  each 
value  of  U.  This  Is  similar  to  nested  FORTRAN  DO  loops.  At  the  comple¬ 
tion  of  each  Innermost  variation,  the  electric  field  components  will  be 
printed  and.  If  specified,  plotted  using  the  scales  specified. 

The  dependent  axis  will  be  the  magnitude  of  the  electric  field 
In  volts/meter  and  the  independent  axis  will  be  one  of  the  geometric 
variables.  Note  that  to  specify  a  polar  plot  with  either  R,  X,  Y,  or  Z 
as  the  most  rapidly  varying  coordinate  is  meaningless  and  will  result  In 
an  error.  However,  requesting  a  linear  or  log  Independent  axis  for  an 
angular  coordinate  Is  not  meaningless  and  will  be  plotted.  A  combination 
R,  T,  P  will  Imply  a  spherical  coordinate  system  while  R,  T,  Z  Implies  a 
cylindrical  system  and  X,  Y,  Z  implies  a  Cartesian  system.  These  are  the 
only  combinations  allowed  and  the  meanings  of  the  primary  coordinate 
Identified  (R,  T,  P,  X,  Y,  Z)  are  given  in  the  table. 

KEYWORD  DEFINITIONS  FOR  COORDINATE  SYSTEMS 


SPHERICAL 

CYLINDRICAL 

CARTESIAN 

R 

Radius  of  Point 
from  Origin 

Radius  in  XY 

Plane 

T 

Angle  measured 
from  Z  axis 
toward  X-Y  plane 

Angle  measured 
from  X  axis 
counterclockwise 

P 

Angle  measured 
from  X  axis 
counterclockwise 

X 

X  Coordinate 

Y 

Y  Coordinate 

Z 

Z  Coordinate 

Z  Coordinate 
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EFIELD  (Continued) 


The  default  values  for  Ul,  VI,  Wl,  DU,  DV,  and  DW  are  all  zero. 
The  coordinate  specified  will  take  all  values  from  U1  to  U2  In  steps  of 
DU.  U2,  V2,  and  W 2  are  defaulted  to  Ul,  VI,  and  Wl.  The  only  case  In 
which  a  geometric  parameter  may  be  completely  defaulted  Is  the  R  Item  in 
the  spherical  coordinate  system.  Otherwise,  the  Initial  value  or  the 
final  value  must  be  given  for  each  of  the  coordinate  axes. 

Example: 

FLD  »  EFIELD  (CUR)  LOGPLR  R1  =  36.0  T1  =  0.0,  DT  =  10.0, 

T2  *  90.0,  PI  =  0.,  DP  *  10.0, 

P2  =  360. 

This  will  cause  the  computation  of  the  near  field  at  a  distance 
of  36  meters  In  the  upper  hemisphere.  Log  polar  plots  of  Er,  E0,  and  E^ 
will  be  made  for  theta  angles  of  0°  to  90°  In  Increments  of  10°  with  37 
points  for  each  component  per  plot. 

The  plot  output  associated  with  the  electric  field  output  is 
unlabeled  In  version  3  of  GEMACS.  The  plots  are  Intended  to  be  used  in  a 
qualitative  manner,  and  the  data  are  listed  prior  to  the  plot.  The 
abscissa  (X)  axis  Is  printed  across  the  page  and  the  ordinate  (Y)  is 
printed  down  the  page.  In  the  event  of  a  polar  plot,  the  origin  is  at 
the  center  of  the  display  region  while  for  nonpolar  plots,  the  origin  is 
as  Indicated  by  the  axis  data  listed. 

All  plots  use  the  most  rapidly  varying  geometric  parameter  as 
the  Independent  variable  (angle  for  polar  plots,  abscissa  for  nonpolar 
plots).  The  dependent  variable  Is  the  ratio  of  the  magnitude  of  the 
electric  field  at  the  observation  point  to  the  maximum  electric  field 
computed  for  an  observation  points.  The  LOG  specification  results  in 
the  value  being  modified  to  20  log-jg  of  the  ratio  (Field  Strength/Maximum 
Field  Strength).  The  dynamic  range  of  the  plots  is  100  dB. 

When  the  Independent  variable  Is  an  angular  coordinate  for  a 
polar  plot,  the  location  of  the  reference  depends  on  the  coordinate  and 
the  coordinate  system  In  use.  For  cylindrical  coordinates,  the  0  angle  is 


EFIELD  (Continued) 


Measured  positive  counterclockwise  from  the  positive  X-axis.  The  same 
convention  holds  for  the  t  angle  In  spherical  coordinates.  This  Is 
Illustrated  In  figure  a.  For  the  spherical  0  angle,  the  measure  Is 
positive  counterclockwise  from  the  plot  X  axis.  In  this  case,  the  plane 
of  the  plot  Is  the  plane  containing  the  vectors  r  and  z,  where  r  is  in 
the  direction  of  the  observation  point  and  z  Is  parallel  to  the 
Cartesian  Z-axIs.  This  Is  Illustrated  In  figure  b. 

Angular  coordinates  for  nonpolar  plots  are  treated  the  same  as 
any  other  Independent  variable  and  plotted  as  the  abscissa. 


End  of  Commands 


END 

The  END  card  designates  the  end  of  a  command  Input  deck.  Any 
text  may  be  on  the  same  card  as  long  as  It  Is  separated  from  the  END  by 
at  least  one  blank  or  comma  (i.e.,  no  comment  delimiter  Is  needed).  In 
addition,  the  END  must  be  the  first  field  encountered. 

Examples? 

END 

END  OF  COMMANDS 

END  OF  FCP747  ANALYSIS 

When  an  END  card  Is  encountered,  execution  Is  terminated,  all 
files  are  closed,  and  an  end-of -module  checkpoint  Is  written. 


ESRC 


Electric  Field  Excitation 


SOS  =  ESRC  t (DS) ]  (,FRQ  =  DSN]  SW  =  DSN,  DSN 
[ ,R  =  DSN]  [, THETA  =  DSN]  [ .PHI  =  DSN] 

[  ,ECC  =  DSN] 

This  command  generates  or  modifies  the  excitation  specified  in 
SDS  by  driving  the  structure  identified  by  DS  with  an  incident  electric 
field.  The  default  structure  identifier  GEODAT  is  supplied  if  omitted. 
The  frequency  in  MHz  is  specified  in  the  FRQ  item;  and,  if  omitted,  the 
value  used  is  the  last  value  specified  in  a  FRQ  item  of  any  command.  If 
the  frequency  has  changed  since  the  last  excitation  (either  ESRC  or 
VSRC),  the  data  are  reinitialized  to  zero  before  the  excitation  is 
computed.  Once  computed,  the  excitation  is  superpositioned  with  previous 
excitation  data.  The  angular  coordinates  of  the  source  are  illustrated 
in  the  figure.  THETA  is  measured  in  degrees  from  the  Cartesian  Z  axis 
and  PHI  is  measured  in  degrees  counterclockwise  from  the  Cartesian  X 
axis.  The  default  values  are  THETA  =  90.,  PHI  =0..  If  the  radial  loca¬ 
tion  specified  in  the  R  item  is  positive,  a  spherical  incident  wave  from 
a  source  located  at  R,  THETA,  and  P  I  will  be  generated.  If  the  radial 
location  is  omitted  or  i.egative,  a  plane  wave  incident  from  THETA  and  PHI 
will  be  generated.  The  default  value  is  R  =  -1.  The  vector  components 
of  the  source  field  are  specified  as  the  values  in  the  SW  item.  The 
first  value  is  the  component  of  the  field  in  the  spherical  9  direction, 
and  the  second  value  is  the  component  in  the  spherical  $  direction.  If 
the  item  is  SW  =  -1.,  0.  for  example,  it  corresponds  to  a  vertically 
polarized  electric  field  with  an  intensity  of  1  volt/m  if  9  =  90°.  If  an 
elliptically  polarized  incident  field  is  to  be  generated,  the  ratio  of 
the  minor  to  major  axis  is  specified  in  the  ECC  field.  The  default  value 
is  ECC  =  0,  indicating  a  linearly  polarized  wave.  Left  or  right 
polarization  is  denoted  by  the  sign  of  this  parameter.  Left  polarization 
is  positive  while  right  polarization  is  negative. 
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ESRC  (Continued) 


etotal  =  ep  +  Ec 


Field  Excitation  Geometry  and  Definition 

When  a  ground  plane  has  been  specified  by  a  previous  ZGEN 
command,  the  total  excitation  is  the  vector  sum  of  the  source  field  and 
the  reflected  field.  The  reflected  field  is  calculated  using  the  reflec¬ 
tion  coefficient  method  discussed  in  the  GEMACS  Engineering  Manual. 
Since  this  is  applied  only  for  a  ground  plane  previously  specified,  the 
user  must  ensure  that  the  ZGEN  or  DMP  commands  precede  the  ESRC  commands 
in  the  execution  stream. 

Example: 

VANT  =  ESRC  SW  =  0.,1.,  THETA  =  45.,  ECC  =  1 

This  will  drive  the  default  structure  (GEODAT)  with  a  left-hand 
circularly  polarized  plane  wave  incident  from  45°  in  the  XZ  plane  of  the 
structure. 
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ESRC  (Concluded) 


When  ESRC  is  used  to  specify  a  source  for  a  problem  with  only 
GTD  geometry  objects,  a  data  set  is  not  generated.  Rather  the  SDS  is 
used  later  by  an  EFIELD  command  to  specify  the  source(s)  which  generated 
the  incident  field.  As  in  the  MOM  and  hybrid  cases,  more  than  one  ESRC 
command  can  be  used.  However,  when  incident  fields  are  requested  through 
an  EFIELD  command,  the  only  ESRC  sources  which  will  contribute  to  the 
incident  field  are  those  on  cards  prior  to  the  EFIELD  command  but  after  a 
LOOP/LABEL  pair,  geometry  change,  or  frequency  change.  In  the  example, 

FRQ  =  500. 

1  SRC1  =  ESRC  (  )... 

FRQ  =  1000. 

2  SRC1  =  ESRC  (  )... 

3  SRC1  =  ESRC  (  )... 

FLD  =  EFIELD  (SRC1 ) _ 

only  the  last  two  ESRC  commands  generating  SRC1  would  contribute  to  FLD. 


GMOATA 


Generate  a  Structure  Geometry 
GHOATA  [=  S]  [,LU  =  N] 

This  command  calls  the  geometry  processor  to  read  the  geometry 
data  cards  which  follow  the  command  stream  END  card.  Note  that  geometry 
data  do  not  get  read  until  after  all  commands  are  read.  The  default 
symbol  name  of  the  geometry  data  set  S  Is  GEODAT.  The  processor  will 
read  the  data  from  the  FORTRAN  logical  unit  specified  In  the  LU  =  N  item. 
Again,  the  user  must  assure  that  N  Is  a  valid  logical  unit  and  the  file 
must  be  In  card  Image  format.  The  default  Is  the  user's  computer  system 
card  Input  file  and  Is  set  as  variable  LUTASK  In  BLOCK  DATA. 

Note  that  If  the  default  name  Is  used.  It  must  not  be  explic¬ 
itly  referenced  In  a  subsequent  command  (l.e..  It  must  be  defaulted  on 
all  commands  for  the  remainder  of  the  command  stream).  The  only 
exception  Is  when  the  PLOT  command  Is  used. 

Examples: 

GHDATA 

GMDATA  =  FCP747 
GMDATA  =  HUT,  LU  =  11 
ZGEN  FRQ=295.0  ZMATRX=MATZIJ 
ZGEN  GMDATA=FCP747  ZMATRX=FCPZIJ 
ZGEN  GMDATA=HUT  ZMATRX=HUTZIJ 

These  examples  result  In  three  geometry  data  sets  being 
generated  and  Identified  as  GEODAT,  FCP747,  and  HUT,  respectively.  These 
data  sets  are  then  used  to  generate  Interaction  matrices.  The  first  ZGEN 
command  uses  the  default  geometry  data  set  and  generates  an  interaction 
matrix  with  the  name  MATZIJ.  The  remaining  two  ZGEN  commands  use  the 
geometry  data  specified  and  generate  Interaction  matrices  with  the  names 
Identified  In  the  respective  ZMATRX  parameter  field. 


Command  Repetition 


LOOP  {J}  N 

label! 


These  commands  cause  the  commands  contained  between  them  to  be 
executed  N  times.  Loops  may  be  nested  to  any  level  as  long  as  the  total 
number  of  LOOP  commands  does  not  exceed  10.  When  using  nested  loops,  the 
loops  must  be  terminated  from  the  innermost  loop  to  the  outermost  loop. 
The  {2}  field  may  be  an  up  to  six  character  alphanumeric  entry  or  an 
integer. 


$  EXECUTE  SUBSEQUENT  COMMANDS  5  TIMES 


LABEL  1 

LOOP  L00P1 ,  5  $  EXECUTE  SUBSEQUENT  COMMANDS  5  TIMES 

LOOP  L00P2,  10  $  EXECUTE  SUBSEQUENT  COMMANDS  10  TIMES 

$  COMMANDS  EXECUTED  50  TIMES 

LABEL  L00P2  $  INNER  LOOP  TERMINATOR 

$  COMMANDS  EXECUTED  5  TIMES 

LABEL  L00P1  $  OUTER  LOOP  TERMINATOR 
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Matrix  Decomposition 
SDS  *  LUD  (DS) 


This  command  results  in  the  decomposition  by  rows  of  the  matrix 
DS  Into  lower  and  upper  triangular  matrices  identified  by  SDS.  In  most 
cases  the  matrix  DS  will  not  reside  in  main  memory,  and  subsequently, 

there  will  be  two  triangular  matrices  generated.  They  will  be  internally 

Identified  by  symbols  in  which  the  rightmost  three  characters  of  the 

symbol  specified  for  SDS  are  replaced  with  LWR  and  UPR  for  the  lower  and 
upper  triangular  matrices  respectively.  If  SDS  has  three  or  less 

characters,  then  the  matrices  are  simply  LWR  and  UPR.  This  will  be 
transparent  to  the  user  since  any  reference  to  SDS  will  result  In  the 
retrieval  of  the  lower  and  upper  triangular  matrices  as  necessary  for  the 
operation.  However,  the  user  may  reference  the  data  symbolically  by 
using  the  original  symbol  with  LWR  and  UPR  replacing  the  last  three 

characters  of  the  name. 

No  pivoting  will  take  place  during  decomposition.  This  is  due 
to  the  fact  that,  for  EM  problems,  pivoting  Is  usually  of  little  value 
and  can  be  quite  time  consuming  for  matrices  not  stored  In  main  memory. 
Pivoting  may  be  added  In  a  subsequent  version  If  a  need  Is  demonstrated. 
Examples: 

LUDZIJ  =  LUD  (ZIJMAT) 

ZIJMAT  =  LUD  (ZIJMAT) 

IM  =  LUD  (ZIJMAT) 

PRINT  LUDLWR,  LUDUPR 
PRINT  ZIJLWR,  ZIJUPR 
PRINT  LWR,  UPR 

Note:  When  the  matrix  OS  does  not  reside  In  main  memory.  It  will 
still  exist  on  the  original  file.  The  LWR  and  UPR  matrices  will  reside 
on  files  other  than  this  original.  Therefore,  purging  DS  will  not  affect 
the  LWR  and  UPR  matrices.  Also,  the  user  will  retain  DS  if  the  LWR 
and/or  UPR  matrices  are  purged. 


V.v*> 


PLOT 


Plot  Geometry  Command 
PLOT  DS 

This  command  will  generate  a  printer  plot  of  the  geometry  In 
data  set  OS.  Three  projections  will  be  made:  X-Y,  X-Z,  and  Y-Z.  These 
provide  a  visual  summary  of  the  location  of  each  wire  segment  and  area 
patch. 

Examples: 

1.  GMOATA 

PLOT  GEOOAT 


2.  GMDATA=HUT 
PLOT  HUT 

In  the  first  example  the  geometry  was  allowed  to  be  given  the 
default  name  GEODAT.  In  order  to  obtain  a  plot  of  the  data,  the  default 
name  for  the  geometry  data  set  had  to  be  explicitly  referenced.  In  the 
second  example  the  geometry  was  named  HUT.  The  data  In  this  geometry  set 
will  be  plotted  due  to  the  presence  of  the  second  PLOT  command. 
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Print  Data  Command 


PRINT  SOS,  SDS, . . .  SDS 

This  command  allows  the  user  to  obtain  the  entire  contents  of 
each  data  set  specified.  If  the  data  are  complex,  the  output  will 
contain  both  the  real  and  imaginary  components  as  well  as  the  magnitude 
and  phase.  For  complex  data  there  will  be  2  elements  per  line;  and  for 
real  data,  there  will  be  10  elements  per  line.  Printing  out  large 
complex  arrays  can  consume  a  fair  amount  of  paper,  and  the  user  is 
advised  to  use  the  WRITE  command  when  the  entire  contents  are  not 
required. 

Example: 

PRINT  VOLT,  CURRNT 

This  command  will  cause  the  data  associated  with  symbols  VOLT 
and  CURRNT  to  be  printed  sequentially. 


PURGE 


Purge  Oata  Command 
PURGE  SOS,  SDS,...  SDS 

This  command  will  cause  the  internal  main  storage  or  FORTRAN 
logical  unit  associated  with  the  specified  symbols  to  be  made  available 
for  other  use.  The  data  are  not  retrievable  after  a  PURGE  command.  Due 
to  the  limited  FORTRAN  logical  units,  it  is  recommended  that  symbols  be 
purged  when  no  longer  required.  This  will  also  make  checkpoint  files 
shorter.  Purged  symbols  may  be  referenced  if  the  actual  data  are  not 
required.  This  could  occur,  for  example,  after  a  matrix  has  been  decom¬ 
posed  and  all  that  is  required  are  the  lower  and  upper  triangular 
matrices. 

Example: 

GMDATA 

ZGEN  FRQ  *  123.0  ZMATRX  =  ZIJMAT 
ZIJMAT  =  LUO  (ZIJMAT) 

PURGE  ZIJMAT 

VOLT  =  VSRC,  V  =  l.,0.  SEGS  =  51-60 

BACSUB  ZIJMAT  *  CUR  =  VOLT 

PRINT  CUR,  VOLT 

END  OF  COMMANDS 

{GEOMETRY  DATA} 

END  OF  DATA 

The  PURGE  command  releases  the  file  ZIJMAT  on  which  the  inter¬ 
action  matrix  is  stored.  However,  the  lower  and  upper  triangular 
matrices  are  unaffected  since  they  are  stored  in  files  labeled  ZIJUPR  and 
ZIJLWR.  Thus  the  BACSUB  command  can  be  executed,  and  a  PRINT  {or  WRITE) 
command  would  output  the  data  in  the  decomposed  matrices. 
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RSTART 


RSTART 


Checkpoint  Restart 

r* 


i 

[INPUT 

Igtd 

r 

[LU  =  NJ 

.MODULE  =  ] 

MOM 

* 

[output] 

, CPNUM  =  N 


This  command  Is  used  to  restart  a  job  from  checkpoint.  The 
checkpoint  file  Is  to  reside  on  the  FORTRAN  logical  unit  specified  in  the 
LU  *  N  Item.  The  default  Item  Is  LU  =  7,  which  is  the  same  as  the  CHKPNT 
default  logical  unit.  If  an  alternate  logical  unit  is  specified,  the 
user  is  responsible  for  assuring  that  GEMACS  can  access  the  unit.  The 
MODULE  Item  Is  used  to  assure  that  execution  Is  restarted  with  the 
correct  module  (INPUT,  GTD,  MOM,  OUTPUT).  The  default  item  is 
MODULE  =  INPUT.  The  checkpoint  number  to  be  recovered  is  specified  in 
the  CPNUM  Item  and  is  not  defaulted.  The  value  N  specifies  the  integer 
number  of  the  checkpoint  to  be  recovered.  This  permits  different  opera¬ 
tions  to  be  performed  on  the  data  without  requirements  to  regenerate  all 
the  data.  For  Instance,  If  a  checkpoint  was  written  after  the  inter¬ 
action  matrix  was  generated,  the  user  could  restart  at  that  point  and  use 
a  different  solution  procedure  than  on  a  previous  analysis.  See  the 
WIPOUT  command  for  further  discussion. 

NOTE:  RSTART  should  be  the  first  conmand  encountered  for  the 
restart  run.  If  it  is  not,  all  previous  commands  in  the 
restart  input  stream  are  Ignored. 

Example: 

RSTART  MODULE=GTD,  CPNUM  =  5 

If  a  restart  is  executed  from  the  same  logical  unit  as  the 
checkpoint  was  written  to,  the  checkpoint  file  may  be  overwritten  on  sub¬ 
sequent  checkpoints.  If  the  user  wishes  to  maintain  the  integrity  of  the 
original  checkpoint  file,  the  restart  must  take  place  from  a  different 
logical  unit  (LU),  and  thus  cannot  be  defaulted. 


Data  Initialization  and  Modification 


SET  SDS  =  N  [ ,N]  [ ,R1  =  HJ  [,R2  =  N]  [ .Cl  =  N] 
t  »C2  =  N] 

This  command  may  be  used  to  initialize  or  change  the  value  of 
data  associated  with  SDS.  If  the  data  are  complex,  then  SDS  =  N, 
N  (corresponding  to  the  real  and  imaginary  components)  is  used.  If  the 
data  are  real,  then  SDS  =  N  is  the  correct  form.  The  items  Rl,  R2,  Cl, 

and  C2  specify  the  row  and  column  limits  of  the  data  to  be  loaded  with 

the  value  specified.  The  default  value  of  Rl  and  Cl  is  1,  while  R2  and 
C2  default  to  Rl  and  Cl  respectively. 

This  command  would  allow  the  user  to  alter  an  excitation  data 
set  if  he  wished  to  force  a  boundary  condition.  If  a  structure  has 
Interior  wires  and  Is  excited  by  an  external  field,  the  field  on  the 

Internal  wires  could  be  reset  to  zero.  Also,  the  Initial  solution  for 

the  BMI  process  could  be  specified,  as  well  as  modifications  made  to  the 
Interaction  matrix. 

Example: 

SET  ZIJMAT  =  0.,0.  Rl  =  10,  C2  =  100 

SET  ZIJMAT  =  0.,0.  R2  =  100,  Cl  =  10 

SET  ZIJMAT  =  l.,0.  Rl  =  10,  Cl  =  10 

SET  VDRV  =  0.,0.,  Rl  =  10 

This  sequence  would  load  zeros  Into  every  element  of  the  tenth 
row  and  tenth  column  of  ZIJMAT  and  then  reset  the  diagonal  element  to 
(1 . ,0) .  The  tenth  element  of  VDRV  would  be  set  to  zero.  This  would  have 
the  effect  of  constraining  the  current  in  the  tenth  segment  to  be  zero 
and  not  allowing  any  interaction  between  the  tenth  segment  and  the  rest 
of  the  structure  except  for  current  continuity  at  junctions. 
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SETINT 


Set  the  Physics  Interactions  for  Matrix  Generation 
SETINT  [KW,  KW,...  KW] 

This  command  selects  which  MOM  and/or  GTD  interactions  are  to 
be  used  to  calculate  interaction,  excitation,  and  field  matrices  and 
whether  or  not  incident  fields  will  be  included  in  the  field  matrix. 
Default  is  the  MOM  interaction.  If  a  SETINT  command  is  not  present  in 
the  input  stream,  only  MOM  interactions  will  be  computed.  Below  is  the 
list  of  Interaction  keywords  and  their  descriptions. 

Wi  Method  of  Moments  interactions 

GTD  All  GTD  interactions  (equivalent  to  PL,  CY,  PC) 

PL  All  plate  interactions  (PR,  PO,  RD,  PDR,  RR) 


PD 

Diffraction  from  plate  edges 

PDR 

Diffraction  from  a  plate  edge;  reflection 
second  plate 

from  a 

PR 

Reflection  from  plates 

RD 

Reflection  from  a  plate;  diffraction  from  a 

plate 

second 

RR 

Reflection  from  one  plate,  then  reflection  from 
another  plate 

All 

cylinder  interactions  (CS,  ER,  ED) 

CS 

Scattering  from  a  cylinder 

ED 

Diffraction  from  end  cap  rim 

ER 

Reflection  from  cylinder  end  cap 

All 

plate-cylinder  interactions  (RC,  rR,  CD,  DC) 

CD 

Scattering  from  a  cylinder,  then  diffraction 
plate  edge 

from  a 

CR 

Scattering  from  a  cylinder,  then  reflection 
plate 

from  a 
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SETINT  (Concluded) 

DC  Diffraction  from  a  plate  edge,  then  scattering  from  a 
cyl inder 

RC  Reflection  from  plate,  then  scattering  from  cylinder 
El  Include  all  excitation  fields  in  Incident  field 

ES  Include  only  fields  scattered  by  GTD  geometry  in  incident  field 

EU  Include  only  fields  not  shadowed  by  GTD  geometry  in  incident 
field 

A  complete  MOM/GTD  hybrid  problem  can  be  worked  with 
SETINT  GTD  MM 

A  MOM  problem  does  not  require  the  SETINT  command,  though  either  of  the 
two  example  commands  could  be  used. 

SETINT  MM 
SETINT 

A  GTD  problem  requires  the  presence  of  the  SETINT  command.  All  keywords 
but  MM  are  valid  in  this  case.  The  example 
SETINT  PL  CY 

would  select  plate  and  cylinder  interactions,  but  not  plate-cylinder 
multiple  interactions. 

The  other  interaction  keywords  provide  the  user  with  a  way  to 
Isolate  scatterers  and  determine  which  types  of  interaction  contribute 
most  to  the  final  results. 
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SOLVE 


Solve  System  of  Simultaneous  Linear  Equations 
SOLVE  QS1*  SOS  -  DS2 

This  command  will  solve  for  the  solution  vector  SOS  using 
lower/upper  triangular  decomposition  on  DS1  and  back  substitution  using 
DS2.  If  DS1  is  already  decomposed,  only  the  back  substitution  wi'l  be 
performed.  SDS  and  0S2  may  be  the  same  symbol. 

Example: 

SOLVE  ZIJMAT  *  CUR  =  E INC 

This  will  cause  execution  of  the  following  equivalent  input: 

ZIJMAT  =  LUO  (ZIJMAT) 

BACSUB  ZIJMAT  *  CUR  =  EINC 
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VSRC 


Voltage  or  Antenna  Excitation 
SOS  =  VSRC  [(OS)]  [ ,FRQ  =  OSN]  ,V  =  DSN,  DSN, 

TAGS 

=  N,  N,  ....  N 

SEGS 

This  command  will  set  up  or  add  to  the  excitation  specified  by 
SDS  on  the  structure  DS.  The  default  structure  identified  is  GEODAT. 
The  voltage  source  is  applied  as  a  delta-gap  electric  field  at  the  mid¬ 
points  of  the  specified  wires  in  the  geometry  data.  That  is,  the 
tangential  electric  field  at  the  midpoints  of  the  wires  specified  is  -V/8, 
where  l  is  the  segment  length  in  meters.  The  frequency  in  MHz  is 
specified  in  the  FRQ  item;  and,  if  the  item  is  omitted,  the  last 
frequency  specified  in  a  FRQ  item  on  any  command  is  used.  If  the  value 
of  the  frequency  has  changed  since  the  last  excitation  (either  VSRC  or 
ESRC),  the  symbol  will  be  reinitialized  to  zero  before  the  new  source 
data  are  computed.  If  the  frequency  is  unchanged,  then  the  source  data 
will  be  added  to  the  existing  data  associated  with  SDS.  This  permits 

superpositioning  of  excitations.  The  real  and  imaginary  components  of 
the  voltage  source  are  specified  by  the  V  item.  Tne  values  for  this 
item  may  have  been  previously  defined  symbolically.  The  segment 

identification  may  have  one  of  two  forms.  If  TAGS  is  specified,  then  all 
wires  which  have  the  tag  numbers  specified  in  the  parameter  list  will  be 
excited.  If  SEGS  is  specified,  only  those  wire  segments  listed  will  be 
excited.  This  list  of  segment  and  tag  numbers  must  be  the  last  entry  of 
the  command  and  may  contain  a  minus  sign  between  successive  entries. 
That  is  N.j,  N-  -  N3,  is  valid  and  will  cause  tags  or  segments 
through  and  including  N^,  and  to  be  excited. 

Example; 

VANT  =  VSRC  ,  FRQ  «  FRQMHZ,  V  --  0.707,  -  0.707,  TAGS  =  1-4 

The  VSRC  command  is  ignored  for  GTD-only  calculations. 
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W I  POUT 


Command  Stream  Modification 

W I POUT  N 

This  command  will  cancel  all  commands  beginning  with  the 
command  identified  by  the  sequence  number  N.  The  commands  as  input  by 
the  user  are  sequenced  in  ascending  order  starting  with  1.  Note  that  the 
command  sequence  number  is  not  necessarily  the  same  as  the  card  number 
containing  the  command  since  comment  cards  and  continuation  cards  may  be 
present  and  are  each  counted.  The  entire  command  sequence  from  N  up  to 
and  including  the  WIPOUT  command  are  eliminated  from  the  execution. 

The  WIPOUT  command  would  normally  be  used  after  a  RSTART 
command  to  change  the  command  sequence.  Additional  commands  may  follow 
the  WIPOUT  command. 

Example: 

RSTART  CPNUM  =  7 

WIPOUT  5 


This  would  eliminate  all  commands  after  the  fourth  command  of 
the  stream  which  generated  the  checkpoint  being  restarted.  Subsequent 
commands  input  after  the  WIPOUT  command  would  be  executed. 

(Note:  LOOP/LABEL  commands  may  be  wiped  out  only  as  paired  sets.) 


WRITE 


Data  Output  Command 

WRITE  DS  [  ,LU  =  N]  [,R1  =  N]  [,R2  =  N]  [,C1  =  N] 

[ ,C2  =  H]  [.FILEIO  =  A] 

Using  this  command,  partial  data  associated  with  the  symbol  DS 
may  be  written  to  the  file  specified.  If  the  logical  unit  item  LU  =  N  Is 
not  specified,  the  system  printer  Is  used  and  the  FILE  ID  item  Is  ignored. 
If  LU  *  N  Is  specified,  the  field  A  In  the  FILEID  item  is  written  as  the 
first  word  on  the  file,  followed  by  the  data  specified.  In  this  case, 
the  output  is  In  FORTRAN  binary  format  and  may  be  used  as  input  to  other 
analyses.  Rl,  R2,  Cl,  and  C2  define  the  row  and  column  limits  of  the 
data.  Default  values  for  Rl  and  Cl  are  1,  while  default  values  for  R2 
and  C2  are  dependent  on  whether  Rl  and/or  Cl  are  specified.  If  Rl  is 
specified,  R2  defaults  to  Rl.  If  not,  then  R2  defaults  to  the  number  of 
rows  in  the  symbol  DS.  The  same  procedure  applies  to  C2. 

Example: 

WRITE  ZIJMAT  Rl  =  10,  Cl  =  1 ,  C2  =  5 

This  prints  the  first  5  elements  in  row  10  of  ZIJMAT 

WRITE  ZIJMAT 

This  writes  the  entire  matrix  ZIJMAT,  and  is  therefore 

equivalent  to  the  results  obtained  by  use  of  the  PRINT  command. 


Interaction  Matrix  Generation 


ZGEN  [.GMDATA  =  DS]  [,FRQ  =  DSN]  , ZMATRX  =  S 
[.ZLOADS  =  DS]  [ ,COND  *  DSN]  [,EPSR  =  DSN] 

This  command  generates  an  Interaction  matrix  using  the  sine  + 
cosine  +  pulse  expansion  function  for  wires,  the  pulse  expansion  function 
for  patches,  and  Impulse  matching  on  the  structure  specified  in  the 
GMDATA  Item.  The  default  structure  Is  GEODAT,  the  default  name  for  the 
GMDATA  command.  The  frequency  In  MHz  is  specified  by  the  value  DSN  of 
the  FRQ  Item.  The  default  Is  the  frequency  given  in  the  last  FRQ  item 
explicitly  stated  In  any  previous  command.  Note  that  the  frequency  may 
be  specified  symbolically  or  numerically.  The  resultant  impedance  matrix 
Is  Identified  by  the  parameter  In  the  ZMATRX  Item.  If  the  structure  is 
to  be  loaded,  the  value  In  the  ZLOADS  Item  must  be  the  symbol  name  of  a 
data  set  generated  by  a  preceding  ZLOADS  command.  The  default  is  no 
loading  or  a  null  ZLOADS  Item.  If  a  ground  plane  is  to  be  used,  the 
conductivity  In  mhos/meter  must  be  specified  In  the  COND  Item.  COND  =  -1 
Implies  a  perfect  ground,  COND  =  0.0  implies  no  ground.  Default  is 
COND  *  0.0.  If  a  nonperfect  ground  Is  specified,  the  relative  dielectric 
constant  er  may  be  specified  in  the  EPSR  item.  The  default  item  is 
EPSR  =  1.  For  a  perfect  or  no  ground  case,  the  contents  of  the  EPSR  item 
are  Ignored.  When  a  ground  is  specified,  it  is  assumed  to  be 
perpendicular  to  the  global  Z  axis. 

Example: 

ZGEN  FRQ  =  FRQMHZ,  ZMATRX  =  ZIJ,  COND  =80.0 

NOTE:  When  the  GTD  interactions  are  specified  on  a  SETINT 
command,  a  ground  plane  may  not  be  specified.  In  this 
case,  both  the  COND  and  EPSR  items  must  be  defaulted. 
The  effect  of  perfectly  conducting  ground  may  be  modeled 
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ZGEN  (Concluded) 


by  locating  one  of  the  GTD  plates  at  z  =  0  and 
specifying  the  PR  keyword  in  a  SETINT  command.  This, 
computationally,  Is  equivalent  to  specifying  COND  =  -1. 


/  / 


ZLOADS 


Structure  Loading 


ZLOAOS  =  SOS  [.GMDATA  =  OS] 


CONO  =  DSN 
ZIMP  =  DSN, 


PRLC 

SRLC 


DSN, 


DSN 

DSN,  DSN 


This  command  allows  a  user  to  place  electrical  loads  on  the 
wire  structure  identified  in  the  GMDATA  item.  The  default  structure 
(GEODAT)  will  be  used  if  the  GMDATA  item  is  omitted.  The  load  informa¬ 
tion  data  will  be  associated  with  the  symbol  SDS. 

The  type  of  loading  is  specified  by  one  of  the  parameters  CONO, 
ZIMP,  PRLC,  or  SRLC.  COND  is  used  to  specify  the  segment  conductivity  in 
mhos/meter.  ZIMP  is  used  to  specify  a  frequency- independent  lumped  load 
resistance  and  reactance  in  ohms.  The  PRLC  and  SRLC  parameters  permit 
loading  wire  segments  with  parallel  (PRLC)  or  series  RLC  (SRLC)  circuits. 
The  values  of  R,  L,  and  C  are  contained  in  the  parameter  list  in  that 
order,  the  latter  two  being  functions  of  frequency.  The  units  of  R,  L, 
and  C  for  the  PRLC  and  SRLC  items  are  ohms,  millihenries,  and  microfarads 
respectively.  The  equivalent  circuits  are  illustrated  below: 


SRLC 
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Z LOADS  (Concluded) 


The  TAGS/SEGS  item  is  used  to  identify  those  wire  segments  to 
be  loaded.  Use  of  the  TAGS  option  results  in  all  segments  which  have  one 
of  the  tags  specified  being  loaded.  Use  of  the  SEGS  option  limits  the 
loading  to  only  those  segments  specified.  The  integer  list  In  the 
TAGS/SEGS  item  may  contain  consecutive  Integers  separated  by  a  hyphen. 
In  this  case,  all  elements  from  the  first  integer  to  the  last  are 
effectively  specified  and  the  corresponding  segments  are  loaded. 

Multiple  ZLOADS  commands  are  permitted.  In  order  to  be 
effective,  the  ZLOADS  command  must  occur  before  the  ZGEN  command. 

Example: 

ZLOADS  »  FCPLOO  GMOATA  =  FCP747,  PRLC  =  5.0,  13.0,  21.0, 

SEGS  -1,5,  7-23,  47 

This  command  causes  the  structure  FCP747  to  have  a  parallel  RLC 
load  applied  to  wires  1,  5,  7,  8,  9,...,  21,  22,  23,  and  47  with  R  *  5 
ohms,  L  -  13  mH,  and  C  -  21  pF. 


Geometry  Input  Language  Processor 


The  function  of  the  GIP  (Geometry  Input  Processor)  is  to  trans¬ 
late  the  user's  Inputs  related  to  structure  geometry  into  a  data  set 
which  may  be  operated  on  by  an  Interaction  matrix  generator  to  provide 
the  interaction  matrix  of  the  structure  under  analysis.  The  GIP  is 
entered  on  encountering  the  GMDATA  command  in  the  task  execution  list. 
It  is  a  re-entrant  processor  in  that  it  may  be  called  several  times, 
either  to  extend  a  previously  generated  geometry  data  set  or  to  create  a 
new  geometry  data  set.  The  attributes  of  the  geometry  data  set  are 
assigned  by  the  GIP.  On  completion  of  the  geometry  processing,  all  data 
are  written  out  to  a  peripheral  storage  device. 

The  user  inputs  to  the  GIP  consist  of  card  images  in  which  the 
first  nonblank  item  is  an  alphabetic  code  designating  the  type  of  data 
contained  on  the  command.  The  data  following  the  type  code  must  be 
separated  by  blanks  or  commas.  Blanks  embedded  within  items  are  not 
allowed.  An  item  is  therefore  defined  as  a  contiguous  group  of 
characters  which,  when  Interpreted,  correspond  to  a  data  input  of  the 
card  being  processed.  Text  following  a  $  is  regarded  as  a  comment  and  is 
Ignored.  There  are  a  maximum  of  256  items  per  command,  not  including 
comment  fields  or  commas  which  are  not  processed.  Continuation  cards  are 
Indicated  In  the  following  ways.  First,  when  a  card  ends  in  a  comma,  the 
next  card  Is  read  as  a  continuation  and  must  have  a  continuation 
character  in  column  1.  Second,  any  card  which  has  a  continuation 
character  In  column  1  Is  assumed  to  be  a  continuation  of  the  previous 
card.  The  only  limit  on  the  number  of  continuation  cards  is  that 
dictated  by  the  limit  of  256  items  per  command.  The  continuation 
character  and  its  definition  are  the  same  as  described  in  section  C. 

The  currently  supported  type  codes  and  their  meanings  are: 

Type  Code  Definition 


Attach 

Combine  Elements 
Connect  Point 
Coordinate  System 


CY 

Cylinder 

0E 

Define  End 

OF 

Define 

EC 

Endcap 

END 

End 

MP 

Multiple  Points 

PA 

Patch 

PL 

Plate 

PSYM 

Plane  Symmetry 

PT 

Point 

RA 

Radii 

RF 

Reflect 

RN 

Renumber 

RSYM 

Rotational  Symmetry 

RX,RY,RZ 

Rotate 

SC 

Scale 

UR 

Wire 

XL 

Translate 

The  command  format  and  use  of  each  of  these  types  are  presented  in  the 
following  section. 

General  guidelines  for  modeling  Include  the  following: 

a.  Segments  must  be  short  compared  to  one  wavelength. 
Lengths  of  0.1  X  should  be  adequate  for  most  purposes.  For  wire  grids 
with  square  mesh,  good  results  have  been  obtained  with  lengths  up  to  0.14 
X.  (See  Volume  II,  GEMACS  Engineering  Manual,  section  G.l.c.). 

b.  Actual  wires  should  be  modeled  with  the  actual  radius. 
Grid  models  should  use  a  wire  radius  about  one-fifth  of  the  segment 
length  In  regions  of  square  mesh.  (See  Volume  II,  GEMACS  Engineering 
Manual,  section  G.l.c  for  more  detailed  comments.) 

c.  Grid  mesh  circumferences  should  not  greatly  exceed  0.5  X. 
Larger  circumferences  lead  to  loop  resonances  and  poor  results. 

d.  Segments  with  lengths  differing  by  more  than  a  factor  of 
two  should  not  be  joined.  Small  angles  (less  than  about  20°)  between 


joined  segments  should  be  avoided.  Unjoined  segments  should  be  separated 
by  a  segment  length  or  more.  The  maximum  number  of  segments  at  a 
junction  Is  limited  to  50.  Segments  are  considered  to  be  joined  when 
their  end  points  are  separated  by  no  more  than  ZERO,  a  variable  set  in  a 
data  statement  in  BLOCK  DATA  BLKDAT.  The  value  of  ZERO  should  be  set  in 
the  code  by  the  user  according  to  his  needs  and  the  limits  of  precision 
imposed  by  his  computer.  It  is  calculated  using  the  following  formula: 

ZERO  =  1/2  *  2'(m_1> 

where  m  is  the  number  of  bits  in  the  mantissa  of  the  computer  system. 

e.  The  maximum  number  of  points,  segments,  defined  elements, 
etc.  that  can  presently  be  input  for  the  geometry  is  discussed  in  section 
D.2  of  this  volume. 

f.  The  renumbering  command  (RN)  permits  the  user  to  specify 
the  geometry  in  the  most  convenient  manner  available  and  to  subsequently 
renumber  the  wire  segments  to  locate  the  near-neighbor  interactions  close 
to  the  diagonal  of  the  interaction  matrix.  The  interaction  terms  between 
the  1th  and  jth  segments  will  be  the  (ij)  and  (ji)  elements  of  the  inter¬ 
action  matrix.  If  the  BMI  bandwidth  chosen  is  b,  then  if  i-j  >  b,  the 
interaction  of  1^h  and  jth  segments  will  not  be  included  in  the  band. 
Too  many  (an  as  yet  unquantified  number)  large  interactions  omitted  from 
the  band  will  cause  the  BMI  solution  technique  to  fail  to  converge. 
Therefore,  numbering  the  problem  such  that  near  neighbors  have  approxi¬ 
mately  equal  segment  numbers  will  cause  large  interactions  to  appear  in 
the  band.  However,  too  many  small  coherent  interactions  outside  the  band 
may  also  cause  the  BMI  solution  technique  to  fail  to  converge.  Note  that 
patches,  plates,  cylinders,  and  cylinder  end  caps  are  to  be  input  after 
all  wire  segments  have  been  defined  and  renumbered.  Therefore,  this 
command  and  its  operation  do  not  apply  to  these  objects. 

g.  Wire  segments  should  be  kept  at  least  0.25  X  from  plate 
edges  for  maximum  accuracy.  Segments  could  be  placed  as  close  as 
0.1  X  if  only  engineering  accuracy  (+  25  percent)  is  required. 


h.  Wires,  but  not  junctions  of  wires,  may  be  connected  to 
patches,  ground,  or  plates,  but  not  to  cylinders  or  end  caps.  A  wire  is 
assumed  connected  if  it  is  within  ZERO  meters  of  a  patch  center,  ground 
plane,  or  within  a  plate  boundary. 

i.  A  plate  Is  assumed  attached  to  a  cylinder  if  the  attach¬ 
ment  edge  is  within  ZERO  meters  of  the  cylinder  surface.  Preferably  the 
plate  should  be  set  slightly  inside  the  cylinder  volume.  Attachment  must 
be  made  with  the  plate  edge  parallel  to  the  cylinder  axis  and  edge 
binormal  perpendicular  to  the  .cylinder  surface  normal  at  the  attachment 
edge. 

j.  Plates  may  not  be  attached  to  end  caps. 

3.  Geometry  Input  Language  Commands 

The  general  form  of  the  commands  available  to  the  user  is: 

TYPE  PI,  P2,  P3,... 

where  TYPE  is  one  of  the  type  codes  listed  in  section  C.2.  PI,  P2,  P3 
are  the  ordered  items  required  or  used  in  the  processing  of  the  command 
specified.  The  items  may  be  separated  by  a  comma  or  a  blank.  Two  of  the 
basic  geometrical  objects  in  the  GIP  are  points  and  wire  segments.  In 
addition,  points  and  wire  segments  may  belong  to  larger  groups.  For 
points,  the  only  larger  group  is  referred  to  as  a  defined  element 
generated  by  a  OF  command,  and  reference  to  such  an  element  will  auto¬ 
matically  reference  all  of  the  points  within  the  element.  Line  segments 
may  also  belong  to  a  group  Identified  by  a  tag  number  in  addition  to  a 
line  segment  number.  The  former  method  (OF)  is  preferred  due  to 
programming  considerations.  Thus,  the  user  may  reference  points,  line 
segments,  or  a  group  of  points  and/or  line  segments.  At  present,  there 
Is  no  way  to  reference  individual  line  segments  unless  they  are  part  of  a 
group;  however,  they  may  be  the  only  object  in  the  group.  Whenever  an 
object  Is  operated  on  to  form  a  new  object,  the  known  attributes  of  the 
source  object  are  automatically  given  to  the  new  object  with  the 
exception  of  group  membership.  This  attribute  will  also  be  transferred 
If  the  group  has  not  been  closed  by  a  DE  (Define  End)  command.  To  the 
GIP,  the  attributes  of  points  are: 

(1)  Point  number 

(2)  Point  location 

(3)  Group  membership. 
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The  attributes  of  a  line  segment  are: 

(1)  Segment  number 

(2)  Segment  tag 

(3)  Location  of  end  points 

(4)  Group  membership 

(5)  Radius  of  wire  segment 

(6)  Segment  connection  data. 

A  third  basic  geometrical  object  is  the  surface  patch.  Within 
GEMACS  it  is  treated  much  like  the  wire  segment.  However,  unlike  the 
wire  segments  which  can  be  grouped  by  tag  number,  each  patch  must  have  a 
tag  number  globally  unique  to  that  patch.  Patches  can  be  part  of  a 
defined  element.  The  only  way  to  reference  an  individual  patch  is  by  tag 
number,  unless  a  defined  element  is  made  up  of  only  one  patch.  In 
version  3  patches  cannot  be  loaded  or  excited  by  a  voltage  source.  They 
can  only  be  treated  as  perfectly  conducting  scatterers.  The  attributes 
of  the  patch  are: 

(1)  Tag  Number 

(2)  Location  of  Center  Point 

(3)  Group  Membership 

(4)  Area 

(5)  Direction  of  Outward  Normal 

(6)  Segment  Connection  Data. 

The  GTD  objects  make  up  the  fourth  type  of  geometry  building 
blocks.  GTD  objects  Implemented  in  version  3  are: 


Code 

Geometry 

Maximum  Number 

PL 

Plates 

14,  each  with  no  more  than  6  corners 

CY 

Cylinders 

1 

EC 

End  caps 

2 

GTD  elements  cannot  be  excited  like  MOM  elements,  and  currents 
cannot  be  found  on  them  directly.  They  are  used  only  to  reflect, 
scatter,  or  diffract  fields. 
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The  attributes  of  the  6T0  elements  are: 

(1)  Plates: 

(a)  User-assigned  number 

(b)  Location  of  corners 

(c)  Group  membership 

(2)  Cylinder: 

(a)  User-assigned  number 

(b)  Major  and  minor  radii 

(c)  Length 

(d)  Coordinate  system 

(e)  Associated  end  cap  numbers 

(3)  End  cap: 

(a)  User-assigned  number 

(b)  Associated  cylinder  number 

(c)  Direction  of  outward  normal. 

In  the  discussion  of  the  commands  that  follows,  it  must  be 
remembered  that  the  Items  are  ordered.  This  is  in  contrast  to  the  items 
of  the  execution  commands  which  were  keyword  indexed  to  achieve  order 
Independence.  Keywords  are  not  used  In  the  GIP  In  order  to  achieve  a 
more  succinct  Input. 

It  should  also  be  noted  that  some  of  the  commands  have  items 
that  may  be  defaulted.  An  example  of  this  is  the  connect  points  (CP) 
command.  In  general,  for  a  command  of  the  form 

TYPE  PI,  P2,  [ ,P3]  [ ,P4]  [,P5] 

In  which  P3,  P4,  and  P5  may  be  defaulted,  it  must  be  kept  in  mind  that 
defaults  can  be  achieved  only  from  right  to  left.  To  specify  P5,  both  P3 
and  P4  must  be  specified.  To  default  P4,  only  P3  must  be  specified, 
while  P5  must  assume  Its  default  value.  When  one  or  more  items  are  to  be 
defaulted,  then  all  items  to  the  right  must  also  be  defaulted. 


:  •  j 

'•  •  'M 
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Attach  Operation 


Parameter  Definition 


AT 


NCS 


Attach  operation  code. 


Point,  tag,  defined  element,  or  plate  to  be  operated 
on. 


An  Integer  for  points,  plates,  or  tag  identifiers; 
the  element  nanre  for  defined  elements. 


Integer  Identifier  of  previously  defined  coordinate 
system  to  which  the  object  or  element  is  to  be 
attached. 


This  command  results  In  a  group  of  points,  segments  and/or 
plates  to  be  translated  and  rotated  into  the  coordinate  system  specified. 
No  new  objects  or  elements  will  be  generated.  This  allows  the  user  to 
input  parts  of  the  geometry  In  local  coordinate  systems  and  then  relocate 
them  to  their  actual  position. 

Examples: 

AT  PT  10  3 

This  would  result  in  the  coordinates  of  point  10  being  changed  to 
those  It  would  have  after  being  relocated  to  the  coordinate  system 
3.  If  point  10  originally  had  x,  y,  z  coordinates  (0,0,0)  and  the 
origin  of  coordinate  system  3  were  located  at  (10,0,0)  with  regard 
to  the  global  system,  then  point  10  would  have  (10,0,0)  as  its  coor¬ 
dinates  after  the  operation. 


AT  (Concluded) 


AT  OF  SPHERE  2 

All  parts  of  SPHERE  would  have  their  coordinates  modified  to  those 
they  would  have  if  they  had  originally  been  defined  in  coordinate 
system  2. 

AT  PL  10  2 

Plate  number  10  would  have  its  corner  coordinates  relocated  to  coor¬ 
dinate  system  2. 


Combine  Elements  Operation 
CE  name!,  name2,  name3,... 


Definition 

Combine  elements  code. 

Element  name  by  which  resultant  group  will  be  known. 

name2,name 3...  Names  of  elements  to  be  combined  into  element  name. 

Note  that  these  elements  must  have  been  generated  by 
a  Define  Element  (OF)  operation  and  will  not  be  avail¬ 
able  under  their  original  names  after  this  operation. 

If  the  user  has  several  elements  that  have  been  defined  by  a  Define 
Element  (DF)  operation,  he  may  combine  them  into  one  element  with  this 
operation.  This  is  useful  when  a  user  has  a  collection  of  generic  shapes 
and  he  needs  to  put  them  together  to  form  another  object.  He  may  then 
combine  all  of  the  elements  under  one  name  for  ease  of  future  reference. 
Example: 

Assume  all  previous  data  are  present. 

AT  DF  BOOM  3 

AT  DF  DISH  1 

CE  SATLIT  CYLNDR  DISH  BOOM 

AT  DF  SATLIT  5 

Coordinate  systems  3  and  1  may  have  been  defined  such  that  the  first 
two  AT  commands  move  elements  BOOM  and  DISH  to  locations  on  element 
CYLNDR.  Then  CYLNDR,  DISH,  and  BOOM  are  combined  into  element 
SATLIT  and  the  resultant  element  positioned  in  coordinate  system 
number  5.  As  can  be  seen,  this  has  eased  the  input  of  geometry 
considerably  when  the  relative  locations  of  defined  elements  are 
known  or  easily  calculated.  If  the  user  wanted  to  define  the 
members  of  BOOM,  DISH,  and  CYLNDR  in  coordinate  system  5,  a  consid¬ 
erable  amount  of  precomputation  would  have  been  required  if  the 
elements'  dimensions  were  on  separate  drawings. 


Parameter 

CE 

namel 


CP 


Parameter 

CP 

N1 

N2 

NSEG 


NT  AG 


NRAD 


Definition 
Connect  points  code. 

Integer  number  of  first  point. 

Integer  number  of  second  point. 

Number  of  segments  to  be  generated  between  N1  and  N2. 
(Default  =  last  NSEG  parameter  from  any  previous 
operation.) 

Tag  number  Identifying  all  segments  between  N1  and 
N2.  (Default  =  last  NTAG  parameter  from  any  previous 
operation.) 

Integer  Index  to  radii  table  entry.  (Default  *  last 
NRAD  entry  from  any  previous  operation.) 


After  the  user  has  defined  the  points  N1  and  N2  (PT  command), 
this  operation  will  connect  NSEG  segments  between  these  points.  Each 
segment  generated  will  have  a  TAG  number  specified  by  NTAG  and  a  radius 
retrieved  from  the  NRAD  entry  of  the  radii  table.  The  default  values  for 
NSEG,  NTAG,  and  NRAD  are  those  values  left  over  from  any  previous  opera¬ 
tion  for  which  they  were  defined.  This  Implies  that  they  must  be  defined 
on  the  first  command  which  requires  them. 

Examples: 

CP  1  2  3  0  5 
or 

CP  1,2, 3,0, 5 

This  would  result  In  generating  three  wire  segments  between 
point  1  and  point  2.  They  would  have  a  TAG  of  0  and  a  radius 
extracted  from  the  fifth  radii  entry. 
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VVO 


CP  (Concluded) 


u:- 


CP  1  4 

Assuming  the  first  example  preceded  this  operation,  three  more 
segments  would  be  generated  between  point  1  and  point  4  with  a 
TAG  of  0  and  radius  extracted  from  the  fifth  radii  entry. 
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cv 


.•V  A 

vy  •;  v 


V  -- 


V*'.% 

y*.  i  .  -  -  ■» 


Parameter 


CS 

NCS 

XC.YC.ZC 


RX,RY,RZ 


Coordinate  System  Specification 


CS  NCS  XC  ,  YC  ,  ZC  ,  RX  ,  RY  ,  RZ 


Definition 


Coordinate  system  code. 

Unique  integer  identifier  for  this  system. 


(x,y,z)  location  of  coordinate  system  origin  with 
respect  to  the  global  coordinate  system  origin. 


Rotation  angles  in  degrees  about  the  x,  y,  or  z  axis 
of  the  global  coordinate  system.  Only  one  angle  may 
be  nonzero.  Negative  angle  specification  results  in 
clockwise  rotation  when  looking  toward  origin  along 
the  axis  of  rotation. 


This  command  permits  the  user  to  specify  additional  coordinate 
systems.  The  NCS  item  on  other  commands  references  the  coordinate  system 
identified  by  this  number.  When  the  NCS  parameter  is  specified  on  other 
commands,  the  transformation  from  or  to  this  coordinate  system  will  be 
made.  This  command  must  precede  all  commands  referencing  this  NCS. 
EXAMPLE: 

CS  3  0.0  2.0  1.0  -30.0  0.0  0.0 


Parameter 


Cylinder  Specification 


CY  n  XR  YR  LEN  [HCSJ 


Definition 

Cylinder  specification  code. 

User-assigned  cylinder  number. 

X-axis  radius  In  the  x-y  plane. 

Y-axis  radius  In  the  x-y  plane. 

Cylinder  length  (z-axis). 

Integer  Identifier  of  previously  defined  coordinate 
system  in  which  the  cylinder  Is  centered.  (Default 
Is  the  last  previous  NCS  entry  on  any  command.) 


This  c  ...miand  allows  the  user  to  generate  a  GTD  elliptical  cylinder 
at  the  origin  of  the  NCS  coordinate  system.  Cylinder  length  Is  from  the 
center  (X  *  Y  *  0)  of  one  end  cap  to  the  center  of  the  other  end  cap,  and 
It  Is  centered  on  the  x-y  plane  of  the  NCS  coordinate  system.  End  cap 
attributes  are  entered  separately  with  the  EC  command.  Version  3  permits 
only  one  cylinder  per  geometry  data  set. 


COORDINATE 
SYSTEM  NCS 


Parameter  Definition 

DE  -  Define  end  code. 

This  operation  ends  or  closes  the  group  identified  in  the 
current  Define  Element  (DF)  operation.  The  group  may  not  be  extended 
except  by  a  CE  operation.  Ail  points,  plates,  patches  and  wire  segments 
generated  since  the  last  DF  operation  belong  to  the  group  identified  by 
the  DF  name.  Points  in  this  group  may  be  referenced  without  regard  to 
group  membership,  however,  the  segment  and  patch  data  may  only  be 
referenced  by  identifying  the  group.  It  is  advisable  not  to  generate 
points  under  a  DF  operation  since  the  storage  available  for  points  is 
more  restricted  than  that  for  segments  and  operations  performed  on  the  DF 
element  will  generate  more  points  which  are  not  usually  required  (see  the 
example  for  the  DF  operation). 


DF 


Define  Element  Operation 
DF  name 

Parameter  Definition 

DF  -  Define  element  code. 

name  -  A  six  character  or  less  identifier  must  begin  with 

alpha  character,  but  can  Include  Integers. 

All  points,  plates,  patches  and  segments  generated  between  this 
command  and  the  next  DE  command  will  belong  to  the  element  "name." 
Nested  DFs  are  allowed.  Note  even  though  points  may  belong  to  an 
element,  they  may  still  be  referenced  Individually. 

Examples: 

DF  BOX 

PT  1  1.0  0.0  0.0 
PT  2  0.0  1.0  0.0 
PT  3  -1.0  0.0  0.0 
PT  4  0.0  -1.0  0.0 
CP  1,  2,  1,  0,  1 
CP  2,  3 
CP  3,  4 
CP  4,  1 
DE 

AT  DF  BOX  1 

Points  1,  2,  3,  and  4,  and  segments  generated  it 
these  points  are  Identified  as  belonging  to  BOX.  All  of  tnest  .j'ncs  -jihJ 
segments  are  then  moved  to  coordinate  system  1.  If  an  ope^tir  ov:- 
performed  on  BOX  which  Involved  generation  of  additional  e’arun::*,  ■-!.*- 
points  belonging  to  BOX  would  also  be  used  as  sources  for  xj..  il  wv.1 
points.  For  this  reason,  use  of  points  In  a  defined  elerntts:  id  Le 

avoided  whenever  possible,  since  the  storage  available  for  poirt?  is  more 
restricted  than  that  for  segments. 

Cylinders  and  end  caps  are  not  allowable  DF  objects 


Parameter 

EC 

n 

+ncyl 


0 

♦ 


Definition 

End  cap  specification  code. 

User-assigned  end  cap  number. 

Cylinder  to  which  end  cap  Is  attached:  +  Indicates 
top  end  cap,  -  Indicates  bottom  end  cap. 

Polar  angle  of  end  cap  outward  normal. 

Azimuthal  angle  of  end  cap  outward  normal. 


This  command  assigns  an  end  cap  to  the  cylinder  ncyl  In  the 
geometry  data  set.  Two  EC  commands  are  required  for  each  CY  command. 
Slanted  end  caps  are  permitted  and  slant  angle  Is  determined  from  the  end 
cap  normal.  They  need  not  be  Identical  for  both  end  caps.  Version  3  of 
GEMACS  limits  the  number  of  end  caps  to  two  and  requires  that  the  end  cap 
normal  be  confined  to  the  (x-z)  plane  ($  =0).  The  end  cap  Is  auto¬ 
matically  located  In  the  cylinder  coordinate  system. 

Z 


End  of  Geometry 


END 


Parameter  Definition 

END  -  End  of  current  geometry  designator. 

This  command  causes  the  GIP  to  stop  reading  input.  It  wii 
then  look  for  wire  junctions,  identify  all  multiple  segments,  print  out 
the  geometry  data  set,  and  write  the  geometry  data  set  to  the  use.'* 
specified  data  set  name. 

Multiple  segments  are  defined  as  those  segments  lying  on  an 
axis  of  rotation  or  in  a  plane  of  reflection.  They  are  identical 
segments  with  the  same  end  points  as  the  generating  segment.  The 
generated  segments  do  not  enter  into  the  interaction  matrix  calculation, 
and  they  are  identified  in  the  segment  data  output  by  a  zero  in  the  NSEG 
column. 


Multiple  Point  Connection  Operation 


MP  NPTS,  HPT,  NP2,  NP3,...,  NPNPTS  [,NSEG  ]  [ ,NTAG]  [ ,HRAO] 


Parameter 


NP1.NP2,, 

NPNPTS 


Definition 

Multiple  point  connect  code. 

Integer  value  of  the  number  of  points  to  be 

connected . 

Integer  identification  of  points  to  be  connected. 
There  must  be  NPTS  of  these  values. 

Number  of  segments  between  each  pair  of  points. 
There  will  be  (NPTS-1)*NSEG  segments  generated. 
(Default  NSEG  =  last  NSEG  value  from  any  previous 
command.) 

Integer  tag  number  Identifying  all  segments 

generated.  (Default  NTAG  *  last  NTAG  value  from  any 
previous  command.) 

Integer  value  of  location  of  wire  radius  in  radii 

table.  (Default  NRAD  =  last  NRAD  value  from  any 

previous  command.) 


This  command  permits  the  user  to  generate  a  set  of  points  to  be 
connected  with  an  equal  number  of  segments  between  each  pair  of  points  as 
listed.  There  Is  no  restriction  on  the  location  of  the  points  and  NPTS 
must  be  greater  than  1. 

Example: 


MP  6  1 


10  5  4  2  0  1 


This  command  would  connect  points  1,  3,  7,  10,  5,  and  4  with  a 
wire  whose  radius  Is  stored  In  the  first  entry  of  the  radii 
table.  There  would  be  two  segments  between  each  pair  of 
points,  and  all  segments  would  have  a  tag  of  zero. 
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PA 


Patch  Specification 


PA 


Parameter 


Definition 


PA  -  Patch  specification  code. 

PT  -  Two-letter  code  Indicating  that  the  patch  is  to  be 

placed  at  a  point,  or  that  points  will  be  used  to 
define  three  corners  of  the  patch. 


N,  N1 ,  N2,  N3 
x,  y,  z 
0 
♦ 


Integer  point  Identifier. 

X,Y,Z  location  of  the  center  point  of  the  >  = 

The  polar  angle  for  the  normal  vector  of  the  v-  ”U. 

The  azimuthal  angle  for  the  normal  vector  of  the 
patch. 


AREA  -  Surface  area  of  patch  in  square  meters  or  user  scaled 

units. 


NTAG  -  Integer  tag  Identifier  of  each  patch,  unique  for  each 

patch. 

HCS  -  Integer  Identifier  of  coordinate  system  {!Vf&ul i  *>- 

previous  NCS  entry  of  any  operation.)  So  ' 


This  command  allows  the  user  to  specify  the  xvy,z  !v  ‘ten  fh 
center  point  of  the  patch,  or  to  place  a  patch  at  a  prevlo  si/  defined 
point.  The  normal  vector  of  the  patch  Is  defined  by  th=  p.;la r  ;ng:s  * 
and  the  azimuthal  angle  $  as  shown  In  the  figure,  or  the  ic-mal  is  den- 
tlfled  by  three  previously  specified  points  at  Its  corners  given  in  a 
counterclockwise  direction  looking  In  the  direction  of  .  e  cut*.ar 
normal  of  the  patch.  The  specification  of  the  patch  ar^i,  s  c  iq  .e  ta 
number  and  an  optional  coordinate  system  complete  the  patu. 
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PA  (Continued) 


z 


Patch  Geometry 


NOTES: 

1)  The  NCS  Item  Is  not  applicable  when  the  PT  Input  format  Is  used 
to  define  the  patch  center.  All  specifications  are  then  given 
In  the  global  coordinate  system. 

2)  A  wire  can  be  connected  between  two  patches.  However,  only  one 
end  of  a  wire  segment  may  be  connected  to  a  patch.  Hence  the 
wire  connecting  the  two  patches  must  consist  of  at  least  two 
wire  segments. 

Examples: 

1)  PA  0.0  0.0  5.0  0.0  0.0  100.0  1 

This  command  would  create  a  surface  patch  parallel  to  and  5 

meters  above  the  x-y  plane  In  the  global  coordinate  system. 

2 

The  patch  has  a  surface  area  of  100  m  ,  and  It  Is  Identified  by 
the  tag  number  1. 


HA  (Concluded) 

2)  PT  4  0.0  5.0  0.0 

PA  PT  4  90.0  90.0  100.0  2 

2 

These  commands  place  a  patch  with  an  area  of  100  m  and 
Identified  by  tag  2  at  point  number  4,  which  has  the  (x,  y ,  z) 
coordinates  of  (0.0,  5.0,  0.0).  The  lormal  vector  of  the  patch 
is  the  y  axis. 


PT 

1 

0.0 

0.0 

0.0 

PT 

2 

-1.0 

0.0 

0.0 

PT 

3 

0.0 

1.0 

0.0 

PT 

4 

1.0 

0.0 

0.0 

PT 

5 

0.0 

-1.0 

0.0 

PA 

PT 

1 

PT 

2  '  3 

PA  PT  1  PT  2  3  4  2.(  1 

2 

These  commands  place  a  2  m  patch  in  the  x-y  plane  centered  at 

t 

the  origin.  The  vertices  of  the  pati  h  are  on  the  x,  y,  -x  and 
-y  axes.  The  normal  vector  of  the  oatch  is  the  z  axis.  The 
same  patch  would  result  if  the  last  command  of  the  inout  stream 
is 

PA  0.0  0.0  0.0  PT  2  3  4  2.0  1 


PL 


Plate  Specification 


PL  n 
or 

PL  n 


N1  N  2  N3  [N4]  [N5]  [N6] 

XI  VI  Z1  X2  Y2  Z2  X3  Y3  Z3 
[X4  Y4  Z4]  [X5  Y5  Z5]  [X6  Y6  26] 


Parameter 

PL 

n 

Nl,  N2 . N6  - 

(X1.Y1.Z1)...  - 

(X6.Y6.Z6) 


Definition 

Plate  specification  code. 

User-assigned  plate  numbe-. 

Point  number  of  plate  corners  (3-6). 

Coordinates  of  plate  corners  (3-6)  in  global  coordi¬ 
nate  system. 


This  command  allows  the  user  to  cafine  a  GTD  plate  by  either 
referencing  previously  defined  point  number1  or  specifying  corner  coor¬ 
dinates.  The  latter  format  Is  preferred  a;  it  does  not  require  addi¬ 
tional  entries  In  the  point  table.  A  plate  may  have  a  minimum  of  three 
and  a  maximum  of  six  corners.  Fourteen  plates  are  allowed  in  one 
geometry  data  set. 

The  plate  in  the  figure  can  be  defined  by  either  of  the 
following:  . 


Triangle  Plate  Input  Ge< metry 
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PL  (Concluded) 


PT 

1 

1.0 

:a) 

1.0 

0.0 

PL  10 

1.0 

i 

1  <j 

O.C 

PT 

2 

3.0 

1.0 

0.0 

* 

3.0 

»  '» 

j 

PT 

3 

2.0 

2.0 

1.0 

* 

2.0 

2/d 

•v: 

PL 

10 

1 

2  3 

The  order  of  entry  of  the  corners  defines  the  r •; •„ ?  w 

the  rlght-hanc  rule.  The  plate  normal  must  point  outworn  from  tne 
geometry  when  diffraction  interactions  hive  been  selects*:.  This  wi ; ' 
define  which  aigle  Is  the  exterior  angle  for  edge  and  we-:g»  «i!f-fra>:tiu.. 
coefficient  ca  culatlons.  For  structures  which  are  not  strict^  ccnv*?.,, 
there  may  be  a  dilemma  if  more  than  one  source  is  used.  in  voe  Mgu/a 
below,  the  unit  normals  are  correctly  oriented  for  the  source  lcc-'t ic-r ■<< 
given. 


<T 

<<T 

((  •  SOURCE 

\ 

SOURCE 

(a) 

(b) 

However,  had  both  sources  been  present  in  the  same  prob  em,  int;.  oed 
answers  would  have  resulted.  The  user  may  overcome  this  problem  by 
working  problems  (a)  and  (b)  separately  ant  superposing  .he  results. 

A  wire  may  be  connected  to  a  pi  ite.  The  connection  is  deter¬ 
mined  by  examining  whether  one  end  of  the  wire  lies  vithin  the  doma:n  •/ 
the  plate.  Connecting  a  single  wire  segment  to  t  o  plates  1r,  not 
allowed,  nor  Is  connecting  a  wire  to  a  p'ate  at  lass  than  a  20°  angle 
from  the  plate.  A  wire  connecting  two  p  ates  must  consist  or'  at  least 
two  segments. 


PS  YM 


Plane  Symmetry  Specification 
PSYM  N 


Parameter  Definition 

PSYM  -  Plane  symmetry  specif icetion  code. 

N  -  Number  of  planes  of  symmetry. 

This  command  allows  the  user  to  specify  that  the  geometry  is 
made  up  of  identical  units  that  have  been  generated  by  use  of  the  reflect 
command  (RF).  The  presence  of  this  command  causes  GEMACS  to  use  the 
symmetry  of  the  geometry  to  more  quickly  cevelop  the  interaction  matrix 
and  to  solve  the  resulting  matrix  equation. 

Example:  (See  figure) 

RA  0.025 

PT  1  0.0  1.0  0.0 

PT  2  0.0  0.0  1.0 

OF  FIGURE 
CP  1  2  30  1  1 
RF  DF  FIGURE  Y  1 
RF  DF  FIGURE  Z  2 
DE 

PSYM  2 
END 

The  first  three  commands  establish  a  radius  for  the  segments 
ana  the  locations  of  two  points.  The  n?xt  five  commands  define  30 
segments  in  the  first  quadrant  of  the  Y-Z  plane  and  then  reflect  them 
across  the  X-Z  plane  followed  by  a  reflection  of  those  60  segments  ac-oss 
the  X-Y  plane.  The  PSYM  command  specifies  that  there  are  two  planes  of 
symmetry  in  the  geometry  (X-Z  and  X-Y). 


PT 


Parameter 

Definition 

PT 

- 

Point  specification  code. 

NPT 

- 

Integer  identifier  of  this  point.  This  most  be  a 
unique  number. 

X.Y.Z 

- 

(X,Y,Z)  location  of  point  with  regards  to  NCS. 

NCS 

- 

Integer  identif ication  of  coordinate  system  for 
(x,y,z).  (Default  NCS  =  last  defined  NCS  entry  on 
any  command.) 

This  command  is  used  to  specify  points.  The  identifier  NPT 
must  be  globally  unique  (not  simply  unique  to  coordinate  system  NCS). 
The  point  (x,y,z)  will  be  transformed  from  coordinate  system  NCS  into  the 
global  system  before  it  is  stored,  if  NCS  is  non-zero. 


10.0  0.0  0.0  0.0 

1.  10 


This  defines  point  number  2  as  being  at  (l.,l.,l.)  in  coordi¬ 
nate  system  10.  The  point  will  be  transformed  to  the  global 
coordinate  system  before  being  stored,  and  its  global  coordi¬ 
nates  in  this  case  will  be  (1.0,  1.0,  11.0). 
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RA 


Radii  Specification- 


RA  Rl,  R2,  R3. 


Parameter 


R1,...Rn 


Definition 

Radii  table  entry  code* 

Floating  point  values  of  • ; . . s.  ai  • 

radii  table  will  be  loaded  poi'.-.i  y  \<  with 
values. 


Instead  of  appending  the  w.r>. 

command  that  generates  a  wire  segment,  the  user  simp's  <  r.tor  to  tr 
In  the  radii  table.  These  entries  are  looses  uo.^..tSu  i> 
items.  Currently  n  <  10. 


Examples: 

RA  .001  .125  .0067 

RA  .025  1.00  .0003 

Load  radii  entries  1  through  6  with  the  olv 


A1  [A2]  [A3]  [INCTA6]  [NCS] 


Parameter  Definition 

RF  -  Reflection  operation  code. 

!PT  ) 

TG  (  -  Two-letter  code  to  Indicate  point,  tag,  defined 

DF  {  element,  or  plate  Is  to  be  operated  on. 

PL  ) 

(  °  ) 

In  I  -  For  PT,  PL  and  TG,  the  integer  identifier;  for  DF, 

)  name  /  the  alphanumeric  name  of  the  defined  element. 


A1  -  Alpha  designation  (x,y,  or  z)  of  axis  along  which 

first  reflection  will  occur. 

A2  -  Alpha  designation  of  axis  along  which  second  reflec¬ 

tion  will  occur,  if  one  is  designated. 

A3  -  Alpha  designation  of  axis  along  which  third  reflec¬ 

tion  will  occur,  if  one  is  designated. 

INCTAG  -  Tag  or  plate  number  increment  parameter.  (Default 

INCTAG  =  0,  wires  only;  must  be  specified  for  patches 
and  plates.) 

NCS  -  Integer  identifier  of  coordinate  system  in  which 

reflection  is  to  take  place.  (Default  NCS  =  last  NCS 
value.) 


This  command  causes  the  symmetry  operation  of  reflecting 
through  the  plane  normal  to  the  axis  specified.  Wire  segments  in  the 
planes  of  reflection  are  allowed.  They  are  identified  as  such  and  will 


RF  (Continued) 


not  enter  into  the  Interaction  matrix  calculation.  However,  patches  In 
the  reflection  plane  are  not  allowed.  Plates  in  the  reflection  plane  are 
allowed  but  the  duplicated  plate  serves  no  purpose.  If  the  reflection 
operation  takes  place  on  an  element  currently  being  defined,  all  segments 
generated  will  be  associated  with  the  element  being  defined.  When 
reflecting  a  patch  or  plate  the  increment  must  be  specified.  The  user 
should  use  care  to  ensure  that  the  new  patch  tags  and  plate  numbers  are 
globally  unique.  In  addition,  the  user  may  not  apply  this  command  to  a 
GTD  cylinder  and  Its  end  caps,  or  to  a  geometry  in  which  a  ground  is 
present. 


Examples: 

RF  OF  WHEEL  X  10 

This  will  cause  all  points,  segments,  patches,  and  plates 
associated  with  WHEEL  to  be  reflected  through  the  y-z  plane  and 
their  tags  to  be  Incremented  by  10.  A  more  complete  geometry 
generating  stream  could  be  as  follows: 

OF  WHEEL 

PT  1  0.  0.  0. 

PT  2  1.  0.  0. 

RZ  PT  2  2  45.0 

RF  OF  WHEEL  X  Y  10 
DE 

This  would  generate  the  points  in  the  positions  shown  in  the 
figure.  Note  that  points  in  the  plane  of  symmetry  are 
regenerated.  This  Is  strictly  for  the  user's  convenience  in 
keeping  up  with  the  point  Identifiers.  NPTS  will  result  In 
2n*NPTS  points  Incremented  sequentially  where  n  Is  the  number 
of  planes  of  reflection  and  NPTS  is  the  number  of  points  being 
reflected.  The  additional  commands: 


RF  (Continued) 

HP  13  1  2  3  1  4  7  1  6  15  1  16  11 

*110  1 

CP  11  10  1  0  1 

CP  3  4  1  0  1 

CP  7  6  1  0  1 

CP  15  16  1  0  1 

would  generate  the  outline  of  a  spoked  wheel.  The  user  could 
have  substituted  5,  9,  or  13  for  point  1,  10  for  2,  8  for  4,  14 
for  6f  and  12  for  16  with  the  same  result. 

2.  OF  WHEEL 

PT  1  0.  0.  0. 

PT  2  1.  0.  0. 

RZ  PT  2  1  45. 

PT  4  0.  1.  0. 

MP  6  1  2  3  1  4  3  1  0  1 

DE 

RF  OF  WHEEL  X  Y  5 

These  commands  would  generate  not  only  Identical  points,  but 
Identical  wire  segments.  These  would  eventually  be  detected 
during  the  search  for  multiple  junctions  and  flagged  as  null 
segments.  The  results  would  be  Identical  to  those  obtained  In 
example  1. 


tissa! 
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Renumber  Operation 


RN  11,  12,  13,...,  -In,  In+2,... 

,  I 

Parameter 

Definition 

RN 

Renumber  operation  code. 

11,. ..I 

Integer  numbers  to  control  the 
segment  numbers. 

resequencing  of  wire 

When  using  the  banded  matrix  iteration  technique.  It  may  be 
Important  to  have  the  wire  segments  numbered  correctly  In  order  for  the 
solution  to  converge.  When  a  RN  operation  is  encountered,  all  wire 
segments  which  have  been  generated  since  the  last  RN  are  renumbered 
according  to  the  sequence  specified  by  the  integers  on  the  RN  command. 
The  resequencing  will  start  with  the  first  segment  generated  since  the 
last  RN  operation.  The  segment  numbers  are  changed  to  correspond  to  the 
Integers  on  the  RN  directive.  When  a  negative  integer  Is  encountered 
(e.g.,  -In),  the  sequence  number  of  the  next  wire  segment  will  be  tl\ 
absolute  value  of  the  negative  Integer,  and  the  sequence  numbers  will  be 
Incremented  by  1  until  the  number  of  segments  identified  by  the  next  item 
has  been  sequenced  (e.g.,  In+2). 

Example: 

Suppose  the  user  has  generated  segments  1  through  10  and  wishes 
to  renumber  the  segments: 

RN  2  4  -5  ,  5  ,  1  ,  3  ,  10 


Old  Segment  Numbers 
New  Segment  Numbers 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

2 

4 

* 

± 

LL 

LI 

9 

LI 

LL 

10 

Neither  GTD  geometries  nor  patches  are  affected  by  the  RN 
command . 


Rotational  Symmetry  Specification 


RSYM  N 


Parameter  Definition 

RSYM  -  Rotational  symmetry  specification  code. 

N  -  Number  of  symmetrical  cells. 

cnr.,-jmd  allow?  tr.e  user  rn  snecifv  tiat  the  geometry  is 
itW:  :-0  or  'i.'entv  t.  Tlr.  < hat  hr  -t  b.  ».n  o?  use  of  the  rota- 

t ivi  («X,  rY,  '"-,c  vV;  >- •  -  •  .  ..  causes  GEMACS  t. 

js.-  .  *  .vraaeirv  of  the  qeoBiet? j  to  mor-  qutkiy  I  up  the  interaction 

mat;1.'  mo  to  so«ve  ■  tr  n*s  iting  ^atrtx  oquatioo. 

•iRanpJe:  (See  nqure) 

RA  0.025 

PT  1  C.O  1.0  0.0 

PT  2  0.0  0,0  1.C 

OF  FIGURE 

CP  1  2  30  1  1 

RX  OF  FIGURE  3  90.0  1 

DE 

RSYM  4 
END 

The  first  three  commands  establish  a  radius  for  the  wire 
segments  and  the  locations  of  two  points.  The  next  four  commands  define 
30  segments  In  the  first  quadrant  of  the  Y-Z  plane  and  then  rotate  it 
around  the  X-axis  to  generate  90  more  segments  in  the  Y-Z  plane.  The 
RSYM  command  specifies  that  there  are  four  cell.,  of  symmetry  in  the 
geometry  (one  cell  in  each  quadrant  of  the  Y-Z  plane). 


RX,  RY,  or  RZ 


Rotation  Operation 


I ADD  ANGLE  £ INCTAG]  [NCS] 


Parameter 


Definition 


RX 

RY  J 

—  HO  CaCi  <-• .  •  u/. 

:  o  «  C 

RZ  ( 

PT  ) 

:cdc 

OF  ( 

elem^r:;  •>: 

plat,  t 

PL  ) 

n  ) 

For  PT,  VI 

anu  iu 

n  ( 

OF,  "name" 

is  the 

name  > 

element. 

n 

I  ADD 

ANGLE 

INCTAG 


?•  i  - 

alpha* 


1  W f  i 


Integer  Indicating  the  number  of  additional  objects 
or  elements  to  be  generated. 

Rotation  Increment  about  specified  rotation  axis 
(positive  is  counterclockwise  when  looking  toward 
origin  along  axis  of  rotation). 

Tag  or  plate  number  increment  for  each  additional 
element.  (Default  =  0,  wires  only;  must  be  specified 
for  patches  and  plates.) 


NCS 


Coordinate  system  identifier.  (Default  NCS  =  last 
NCS  specified.) 


The  rotation  operation  can  be  used  to  generate  structures  which 
have  axes  of  revolution.  If  the  objects  to  be  rotated  are  members  of  a 
defined  element,  the  segment  numbers  of  the  original  segments  increment 


RX,  RY,  or  RZ  (Continued) 

by  the  number  of  segments  in  the  original  defined  element  for  each 
additional  element.  If  the  OF  operation  is  still  in  effect,  all 
segments,  patches,  plates  and  points  generated  will  be  members  of  the 
element  being  defined.  When  rotating  a  patch  or  plate  the  increment  must 
be  specified.  The  user  must  ensure  that  the  new  patch  tags  or  plate 
numbers  are  unique. 

This  command  may  not  be  applied  to  a  6TD  cylinder  or  to  its  end 

caps. 


Example: 

RA 

.01 

OF 

CONE 

PT 

Q 

• 

o 

• 

0. 

PT 

2  1.  0. 

10. 

PT 

3  0.  0. 

10. 

RZ 

PT  2  1 

45.  $  ROTATE  PT2  45°  AROUND  Z  AXIS;  GENERATE  PT4 

MP 

3  3  2 

4  1  0  1  $  CONNECT  PTS  3,  2,  4 

CP 

2  1  10 

$  CONNECT  PTS  1  AND  2  WITH  10  SEG  OF  RADIUS  .01 

RZ 

DF  CONE 

7  45.  $  ROTATE  BASIC  ELEMENT  TO  COMPLETE  CONE 

DE 

The  RA  command  establishes  the  radius  of  all  the  segments 
generated.  This  is  followed  by  the  define  element  (OF)  command.  All 
points  and  segments  generated  from  here  to  the  next  define  end  (OE)  card 
will  belong  to  the  structure  called  CONE.  The  next  three  commands  define 
the  basic  points,  the  result  of  which  is  shown  in  figure  a.  Then  point  2 
is  rotated  around  the  Z-axis  to  generate  point  4.  The  next  two  commands 
connect  these  four  points  to  define  the  basic  pattern  of  the  element 
CONE.  This  is  shown  in  figure  b.  This  basic  pattern  is  then  rotated 
around  the  Z-axis  and  regenerated  seven  more  times  to  produce  a  wire- 
gridded  cone.  The  DE  command  then  closes  the  generation  cards  for  CONE. 
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Scale  Parameter 


(  value 

SC 

)  fT  1 

/  IN 

\  CN 

Parameter 


Definition 


SC 


Scale  specification  code. 


Value  Is  the  numeric  value  of  the  scale  factor  In 
meters/unit.  FT,  IN,  CM  are  two-letter  codes  that 
automatically  scale  the  data  from  feet.  Inches,  and 
centimeters  to  meters.  (Default  value  -  1) 


The  GEMACS  code  uses  the  MKS  system  of  units.  Unless  specified 
otherwise,  all  geometry  Input  Is  assumed  to  be  In  meters.  When  an  SC 
command  Is  read,  all  subsequent  dimensions  until  the  next  SC  command  are 
scaled  by  the  Item  on  the  former  SC  command. 

Example; 

SC  FT 

All  subsequent  data  would  be  converted  from  feet  to  meters 
before  being  stored. 


UR  XI,  Yl,  Z1,  X2,  Y2,  Z2,  [,NSEG]  [,NTAG]  [ ,NRAD]  [NCS] 


Parameter  Definition 

MR  -  Wire  Input  designator. 

XT,YT,Z1  -  Coordinates  of  what  will  be  considered  the  negative 

end  of  the  wire  segment. 

X2,Y2,Z2  -  Coordinates  of  what  will  be  considered  the  positive 

end  of  the  wire  segment. 

NSEG  -  Number  of  segments  Into  which  the  wire  Is  to  be 

divided.  (Default  to  last  NSEG  entry  of  any  previous 
operation.) 

NTAG  -  Integer  tag  Identifier  of  each  segment.  (Default  to 

last  NTAG  entry  of  any  previous  operation.) 

NRAD  -  The  location  of  the  radius  of  each  segment  In  the 

radii  table.  (Default  to  last  NRAD  entry  of  any 
previous  operation.) 

NCS  -  Integer  Identifier  of  coordinate  system.  (Default  to 

last  NCS  entry  of  any  previous  operation.) 

This  Is  the  preferred  method  of  generating  wire  segment  data  In 
GEMACS.  It  will  generate  NSEG  segments  between  (XI,  Yl,  Zl)  and  (X2,  Y2, 
Z2)  with  tag  Identifiers  of  NTAG,  of  radius  specified  by  NRAO.  If  NCS  Is 
not  equal  to  zero,  (XI,  Yl,  Zl)  and  (X2,  Y2,  Z2)  will  be  the  end  points 
In  coordinate  system  NCS. 

Example: 

MR  1.63,  2.47,  3.67,  26.4,  16.3,  43.,  10,2,1,1 
This  would  generate  10  segments  with  tag  number  2  and  with  a 
radius  specified  by  the  first  radius  entered  previously  In  the 
radii  table.  The  end  points  would  be  transformed  from 
coordinate  system  1  to  the  global  system  before  storing  segment 
end  points. 


XL 


Translation  Operator 


I  ADO  DX  DY  DZ  [INCTAG]  [NCS] 


Pargetgr 

XL 


IADO 

DX.DY.DZ 

INCTAG 

NCS 


Definition 

Translation  operation  code. 

Two  letter  code  designating  point,  tag,  defined 
element,  or  plate  to  be  translated. 


For  PT,  PL  and  TG,  n  Is  the  Integer  Identifier;  for 
DF,  "name"  Is  the  alphanumeric  name  of  the  defined 
element. 


Integer  Indicating  the  number  of  additional  objects 
or  elements  to  be  generated. 

Incremental  translation  vector.  Each  must  be 
present,  even  If  equal  to  zero. 

Tag  or  plate  number  Increment  for  each  additional 
element.  (Default  »  0,  wires  only;  must  be  present 
for  patches  and  plates.) 

Coordinate  system  In  which  translation  Is  to  take 
place.  (Default  -  last  NCS  value  used  on  any 
previous  command.) 


The  translation  operation  will  generate  I  ADD  additional 
segments,  patches,  points,  plates  and  defined  elements  Identical  to  those 
specified.  Each  will  be  displaced  from  the  original  by  (DX,DY,DZ)  In  the 
coordinate  system  specified.  If  IADO  Is  zero,  the  objects  are  simply 
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XL  (Continued) 

translated  to  a  new  location.  This  last  operation  can  be  performed  by 
defining  a  coordinate  system  at  the  new  location  and  performing  an  AT 
(attach)  operation. 

When  translating  a  patch  or  plate  the  Increment  must  be 
specified.  The  user  must  ensure  that  the  new  patch  tags  or  plate  numbers 
are  globally  unique.  This  command  may  not  be  applied  to  a  GTD  cylinder 
or  Its  end  caps. 


Examples: 

RA 

.01 

PT 

1  0.  0.  0. 

PT 

2  0.  1.  0. 

PT 

3  0.  0.  1. 

DF 

GRIO 

HP 

3  3  12  1 

1  1 

XL 

DF  GRIO  9  0. 

1.  0. 

HR 

0.  10.  0.  0. 

10.  1.  1 

XL 

OF  GRIO  9  0. 

0.  1.  1 

WR 

DE 

0.  0.  10.  0. 

o 

• 

o 

• 

o 

r— 

The  HP  command  would  generate  segments  1  and  2  as  Illustrated 
In  figure  a.  The  first  XL  command  would  generate  segments  3 
through  20,  and  the  WR  command  would  result  In  segment  21  of 
figure  b.  The  second  XL  command  would  generate  segments  22 
through  210,  and  the  UR  command  would  close  the  grid  with 
segments  211  through  220.  This  Is  Illustrated  In  figure  c. 
Future  references  to  GRIO  would  encompass  segments  1  through 
220. 
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Figure  c.  Completed  Wire  Grid 


PA  0.0  0.0  S.O  0.0  0.0  100.0  2 

XL  TG  2  2  0.0  10.0  0.0  1 

This  pair  of  commands  would  generate  patches  located  at  (0.0, 

0.0,  5.0),  (0.0,  10.0,  5.0),  and  (0.0,  20.0,  5.0)  with  tag 

numbers  2,  3,  and  4  respectively.  Each  patch  has  a  surface 

2 

area  of  100  m  and  its  normal  vector  is  In  the  Z  direction. 


D.  COMPUTER  REQUIREMENTS 

GEMACS  Is  written  In  American  Standard  FORTRAN,  X  3.9-1966.  It  Is 
capable  of  executing  with  no  library  subroutines  other  than  those 
required  by  the  ANSI  standard.  The  code  consists  of  four  sequentially 
executable  modules,  whose  computer  memory  requirements  are  given  In 
table  2.  These  values  are  approximate  and  depend  on  the  specific  machine 
and  load  method  utilized. 


TABLE  2.  GEMACS  COMPUTER  MEMORY  REQUIREMENTS 


Module 

Decimal 

Core  Words 

Module 

MOM 

Is  required 
GTD 

for 

Hybrid 

INPUT 

70K 

X 

X 

X 

GTD 

120K 

X 

X 

MOM 

82K 

X 

X 

OUTPUT 

50K 

X 

X 

X 

The  four 

modules  execute 

sequentially,  as 

shown  In 

figure  1. 

However,  should 

the  functions  of 

a  module  not  be  required 

In  problem 

solution,  that  module  may  be  bypassed.  For  Instance,  a  MOM-only  problem 
would  not  require  the  GTO  module,  so  only  the  INPUT,  MOM,  and  possibly 
OUTPUT  modules  would  be  needed. 

The  only  communications  medium  between  the  modules  Is  the  checkpoint 
file.  This  file  contains  the  contents  of  Important  commons  and  the 
contents  of  all  data  files  present  at  the  end  of  a  module's  execution. 
This  allows  two  types  of  runs  to  be  made.  One  may  either  execute  all 
four  modules  In  one  batch  job,  or  one  may  execute  the  modules  In  succes¬ 
sive  jobs,  examining  all  Intermediate  results  before  proceeding  to  the 
next  module.  If  the  latter  Is  the  case,  the  checkpoint  file  should  be  a 
permanent  file. 
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Execution  always  begins  with  the  INPUT  module.  The  INPUT  module 
examines  the  user's  commands  and  checks  for  correctness  and  compatibil¬ 
ity.  The  geometry  Inputs  are  processed  to  generate  the  geometry  data 
sets  required  for  problem  execution. 

The  GTD  module  performs  all  calculations  required  for  geometrical 
theory  of  diffraction  modeling.  Since  It  does  not  generate  any  results 
directly,  the  only  output  from  the  GTO  module  Is  a  series  of  status 
messages  which  Indicate  the  GTO  tasks  being  accomplished. 

The  MOM  module  performs  all  calculations  required  for  method  of 
moments  modeling  and  numerical  solution  of  the  resulting  matrix  problem. 
The  MOM  module  generates  solution  currents  and  Impedances  and  displays 
data  set  contents  via  PRINT  or  WRITE  commands. 

The  OUTPUT  module  is  the  last  module  to  execute.  It  produces  field 
pattern  printouts  and  plots. 

1.  I/O  Requirements 

GEMACS  makes  extensive  use  of  peripheral  file  storage  and  must 
have  several  FORTRAN  logical  units  available.  These  are  listed  in 
table  3  with  the  Internal  variable  name,  the  logical  unit  number ,  and  the 
file  usage  given.  Data  are  stored  starting  on  logical  unit  IOSYNB  up  to 
and  Including  the  logical  unit  number  specified  on  the  NUMFIL  Input.  The 
user  Is  responsible  for  assuring  that  GEMACS  can  access  these  files.  If 
more  files  are  required  than  made  available,  a  fatal  error  will  occur, 
and  an  attempt  will  be  made  to  write  a  checkpoint.  If  a  logical  unit 
(LU)  designator  is  specified  on  the  CHKPNT,  GMOATA,  RSTART,  or  WRITE 
commands.  It  should  be  units  3  or  4,  or  greater  than  that  specified  on 
the  NUMFIL  command.  This  will  ensure  that  the  unit  is  not  In  use  when 
required.  Table  4  lists  the  types  of  data  sets  and  the  number  of  logical 
units  required  for  storage  for  each  type. 

Execution  of  a  PURGE  command  releases  the  logical  unit  used  for 
the  specified  data  set  for  other  use. 
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TABLE  3.  GEMACS  LOGICAL  UNITS 

mmsuammam 

LOGICAL  UNIT 

DESIGNATOR 

NUMBER 

USE 

I0SCR1 

I 

Scratch  File 

I0SCR2 

2 

Scratch  File 

LUTASK 

5 

Card  Image  Input 

LUPRNT 

6 

Formatted  Output 

IOCKPT 

7 

Checkpoint  File  (Binary) 

IOSYMB 

8 

Data  Storage  (Binary) 

LUNIT 

9 

Data  Storage  (Binary) 

LUNIT 

10 

• 

• 

Data  Storage  (Binary) 

LUNIT 

• 

NUMFIL 

Data  Storage  (Binary) 

2.  Internal  Storage  Requirements 

There  are  several  primary  arrays  used  by  GEMACS  for  problem 
execution.  All  arrays  are  stored  in  common  blocks  and  have  internal 
variables  specifying  their  size. 

SEGTBL  Is  used  to  store  the  segment,  patch,  plate,  GTO 
cylinder,  and  end  cap  data.  It  is  dimensioned  as  a  NPRSEG  by  MAXSEG 
array.  The  geometry  data  are  stored  on  a  logical  unit  (LU)  and  as  needed 
are  Input  to  and  output  from  main  memory  in  blocks  of  at  most  MAXSEG 
columns.  Although  MAXSEG  is  presently  set  at  500,  this  does  not  limit 
the  number  of  wire  segments  and  patches  that  are  used  to  model  the  MOM 
geometry.  This  limit  is  20,000  segments  and  patches.  As  stated 
previously,  version  3  of  GEMACS  allows  for  1  cylinder,  2  end  caps  and  14 
plates  which  may  have  at  most  6  corners. 

PTTBLE  Is  used  to  store  point  data  and  Is  dimensioned  as  a 
MAXPTS  by  NPRPT  array.  MAXPTS  is  the  maximum  number  of  points  to  be 
stored,  which  Is  presently  dimensioned  for  100  points. 

IDEFIN  is  used  to  store  the  data  for  defined  elements  and  is 
dimensioned  as  a  MAXOEF  by  NPRDEF  array.  MAX0EF  is  the  maximum  number  of 
defined  elements  allowed,  which  Is  presently  100  elements. 

CVAL  is  used  to  store  coordinate  system  data  and  is  dimensioned 
as  a  MAXCSY  by  6  array.  MAXCSY  is  the  maximum  number  of  coordinate 
systems  allowed,  now  set  to  10.  CVAL  Is  equivalent  to  CX  in  the  CSYSTM 
common  block. 

The  array  TEMP  Is  stored  In  common  TEMP01  and  must  be  of  dimen¬ 
sion  HTEMPS,  presently  set  to  5500.  This  array  is  used  throughout  the 
code  for  Internal  computation  and  intermediate  storage.  At  present,  TEMP 
must  be  capable  of  containing  at  least  three  columns  of  the  interaction 
matrix.  This  Implies  that  NTEMPS  must  be  greater  than  six  times  the 
value  of  MAXSEG  since  the  interaction  matrix  is  complex. 

All  array  parameters  are  stored  in  the  same  common  block  as  the 
array  and  are  preset  In  block  data  subroutine  BLKDAT. 

The  maximum  number  of  plates  allowed  for  GTO  modeling  may  be 
Increased  by  changing  (I)  the  dimensions  of  arrays  associated  with  plate 
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modeling,  (2)  the  variables  used  for  plate  geometry  processing  and  (3) 
the  variables  directing  execution  of  the  GTD  phsylcs.  The  arrays  that 
must  be  changed  are: 


ARRAY  NAME 

XI(14,14,3) 

VXI(3,3,14) 

LSURF(14) 

LSHD(14) 

LIHD(14, 14) 

LDC(14,6) 

FNP(14,6) 

X(14,6,3) 

V(14,6,3) 

VP(14,6,3) 

VH(14,3) 

MEP(14) 

LRFI(14) 

LRFS(14) 

LDRC(14,6) 

LRDC(14,6) 

XX(14,6,3) 


COMMON  BLOCK 
IMAINF 
IMAINF 
SURFAC 
LSHOT 
LSHDT 
LDCBY 
FNANG 
GEOPLA 
GEOPLA 
GEOPLA 
GEOPLA 
GEOPLA 
CLRFI 
CLRFS 
CLDRC 
CLRDC 
PLAIN 


The  plate  geometry  variables  are  located  In  COMMON  /GTDDAT/  and  are 
defined  as  (default  In  parentheses): 

MAXPLT  (14)  maximum  number  of  plates 

The  variables  directing  GTD  physics  are  found  In  subroutine  GTDDRV  (GTD 
module): 

MXPLT— maximum  number  of  GTD  plates 
3.  System  Library  Routines 

No  system  library  routines  are  required;  however,  some  are 
desirable.  The  most  Important  is  a  routine  to  return  the  elapsed 
processing  time  In  minutes.  This  routine  Is  called  from  subroutine 
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TICHEK  and  must  be  available  when  using  the  CHKPNT  command  with  the  CPINC 
item. 

Auxiliary  routines  to  return  the  date  and  time  are  called  by 
subroutine  SYSRTN.  In  the  absence  of  these  routines,  zeros  should  be 
returned  to  the  calling  routine. 

The  file  status  function  routine  LUSTAT  Is  called  after  each 
read  to  detect  an  end  of  file.  If  a  library  function  is  available  to 
determine  this  information.  It  should  be  called  from  this  routine.  If 
none  is  available,  simply  return  a  zero  value  for  the  function. 

4.  Code  Verification  Example 

The  following  example  problem  can  be  used  to  test  that  all  four 
GEMACS  modules  have  been  properly  loaded  on  the  host  computer.  The  INPUT 
module  is  run  first  with  the  commands  and  geometry  data  as  input. 


$ 

$  GEMACS  GTO  TEST  PROBLEM  —  THIELE'S  MONOPOLE  OVER  FINITE  PLATE 
S  DISTANCE  OF  MONOPOLE  FROM  PLATE  EDGE  IS  0.75  METERS 

S  TESTS  STRUCTURE  OF  CODE 
$  MODULE  COMMUNICATIONS 

$  ALL  PLATE  GTD  INTERACTIONS 

S  VALIDATION  OF  SIMPLE  PHYSICS  (NEAR-FIELD) 

$ 

S . 

s 

DEBUG  STATS 
NUMFIL=17 
TIME =4. 

FRq=500. 

SET1NT  PL  MM  El 
GMOATA-THIELE 

ZGEN  GMOATA  -  THIELE  ZMATRX  -  ZGARY 
SRC* VSRC (THIELE)  V*1.,0.  SEGS-1 
SOLVE  ZGARY*I*SRC 
PRINT  I 

FLDOIP*EFIELD(  I )  P1*0.  T1-0.  DT*5.  T2-180. 

END 

RA  .00762 
PT  1  0.  0.  0. 

PT  2  0.  0.  0.25 

CP  1  2  6  1  1 

PT  3  -0.75  -0.75  0. 

PT  4  0.75  -0.75  0. 

PT  5  0.75  0.75  0. 

PT  6  -0.75  0.75  0 

PL  1  3  4  5  6 

END 


The  INPUT  module  will  generate  a  list  of  geometry  elements. 
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WIRE  SEGMENT  DATA 


ENDPOINTS 

CENTERPOINTS 

RAOII  LENGTH  EN01 

1  EN02  NSEG  TAG 

XN 

XP 

YN 

YP 

ZN 

ZP  XC 

YC 

ZC 

0. 

0. 

0. 

0. 

0. 

4.167E-02  0. 

0. 

2.083E-02 

7.620E-03 

4.167E-02 

0 

-2 

1 

0. 

0. 

0. 

0. 

4.167E-02  8.333E-02  0. 

0. 

6.250E-02 

7.620E-03 

4.167E-02 

1 

-3 

2 

0. 

0. 

0. 

0. 

8. 333E-02  0.125  0. 

0. 

0.104 

7.620E-03 

4.167E-02 

2 

-4 

3 

0. 

0. 

0. 

0. 

0.125 

0.167  0. 

0. 

0.146 

7.620E-03 

4.167E-02 

3 

-5 

4 

0. 

0. 

0. 

0. 

0.167 

0.208  0. 

0. 

0.188 

7.620E-03 

4.167E-02 

4 

-6 

5 

0. 

0. 

0. 

0. 

0.208 

0.250  0. 

0. 

0.229 

7.620E-03 

4. 167E-02 

5 

0 

6 

TOTAL  HIRE  SURFACE 

AREA  I.197E-02 

END  1  OF  SEGMENT 

1  IS  CONNECTED 

TO  PLATE  1 

AT  (X.Y.Z)-  0. 

0. 

0. 

PLATE  DATA 

PLATE* 

X 

T 

Z 

X 

Y 

Z 

X 

Y 

Z 

1 

•0.750 

-0.750 

0. 

0.750 

-0.750 

0. 

0.750 

0.750 

0. 

-0.750 

0.750 

0. 

Upon  termination,  the  INPUT  module  writes  an  end-of -module  checkpoint  to 
be  used  In  starting  the  GTD  module. 


CHECKPOINT  NUMBER  1  STARTED  AT  TIME  0.124  ON  LOGICAL  UNIT  7 
COMMON  BLOCK  AOEBUG  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  AMPZIJ  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  ARGCOM  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  CSYSTM  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  DEFDAT  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  FLDCOM  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  GEOOAT  WRITTEN  OUT  TO  ICKFIL 

COMttN  BLOCK  IOFLES  WRITTEN  OUT  TO  ICKFIL 

COMttN  BLOCK  OUNCOM  WRITTEN  OUT  TO  ICKFIL 

COMttN  BLOCK  PARTAB  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCX  SCNPAR  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  SEGMNT  WRITTEN  OUT  TO  ICKFIL 

COMttN  BLOCK  SYHSTR  WRITTEN  OUT  TO  ICKFIL 

COMttN  BLOCK  SYSFIL  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  TEMPO  1  MITTEN  OUT  TO  ICKFIL 

COMttN  BLOCK  6TD0AT  WRITTEN  OUT  TO  ICKFIL 

COMttN  BLOCX  MOOULE  WRITTEN  OUT  TO  ICKFIL 

COMttN  BLOCK  INTMAT  WRITTEN  OUT  TO  ICKFIL 

PERIPHERAL  FILE  B  SYMBOL  THIELE  NUMREC-  1 
CHECKPOINT  COMPLETE  AT  0.131 
ELAPSED  TIME-  0.007 
CURRENT  FILE  LEN6TH-  23421 
STOP  AT  LINE  1 . 


The  GTD  module  begins  execution  from  the  checkpoint.  No  Input 

data  are  used. 
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GEMACS  START-UP  PROCESSOR  CALLED  ON 

8ESIN  START-UP  PROCEOURE 

03/02/83  AT  10.02 

COMMON 

A0E8UG 

READ 

WITH 

1107 

WORDS 

COMMON 

AMPZIJ 

READ 

WITH 

71 

WORDS 

COMMON 

ARGCOM 

REAO 

WITH 

104 

WORDS 

COMMON 

CSYSTM 

READ 

WITH 

72 

WORDS 

COMMON 

OEFOAT 

REAO 

WITH 

506 

WOROS 

COMMON 

FLDCOM 

REAO 

WITH 

7 

WORDS 

COMMON 

GEODAT 

REAO 

WITH 

120 

WORDS 

COMMON 

IOFLES 

REAO 

WITH 

200 

WORDS 

COMMON 

JUNCOM 

REAO 

WITH 

205 

WOROS 

COMMON 

PARTAB 

REAO 

WITH 

3407 

WOROS 

COMMON 

SCNPAR 

REAO 

WITH 

755 

WOROS 

COMMON 

SEGMNT 

REAO 

WITH 

5531 

WOROS 

COMMON 

SYMSTR 

REAO 

WITH 

102 

WORDS 

COMMON 

SYSFIL 

REAO 

WITH 

36 

WOROS 

COMMON 

TEMP01 

REAO 

WITH 

5502 

WOROS 

COMMON 

GTOOAT 

REAO 

WITH 

14 

WOROS 

COMMON 

MODULE 

REAO 

WITH 

13 

WOROS 

COMMON 

INTMAT 

REAO 

WITH 

131 

WOROS 

PERIPHERAL  FILE  8  SYMBOL  THIELE  NUMREC*  1 


START-UP  COMPLETE 

MODULES  PREVIOUSLY  RUM:  INPUT 


The  only  output  from  the  GTO  module  is  a  list  of  status  messages,  as 
shown  below: 

GENACS-GTO  TASK  EXECUTION  STARTEO  ON  03/02/83  AT  10.02 
NUMBER  OF  PERIPHERAL  FILES  AVAILABLE  17 


RUN  TIME  SET  TO 

FREQUENCY  SET  TO 
MAVLENGTH  0.600 


4.00  CPU  MINUTES 

500.  MEGAHERTZ 
METERS 


FILL  IMPEDANCE  MATRIX  ZGARY 
USING  BASIS  FUNCTION  SINCOS 
ON  GEOMETRY  DATA  SET  THIELE 
SHAOOMING  MATRIX  ZGSHO 
FREQUENCY (MEGAHERTZ)  500.00 


INTERACTIONS  SET:  PR  RR 

EXCITE  GEOMETRY  DATA  THIELE 

PO 

RO 

POR 

PL 

MM 

EU 

ES 

EXCITATION  VOLTGE 

EXCITATION  OATA  SRC 

INTERACTIONS  SET:  PR  RR 

PO 

RO 

POR 

PL 

MM 

EU 

ES 

VOLTAGE  EXCITATION  WILL  BE  PERFORMED  IN  MOM  MODULE 

GENERATING  GREENS  FUNCTION  MATRIX  FLOGFM 

MM 

EU 

ES 

INTERACTIONS  SET:  PR  RR 

PO 

RO 

POR 

PL 

GENERATING  INCIOENT  FIELD  MATRIX  FLDOIP 

MM 

EU 

ES 

INTERACTIONS  SET:  PR  RR 

PD 

RO 

POR 

PL 

Finally,  an  end-of-module  checkpoint  Is  taken  at  the  conclusion  of  the  GTD 
module. 

The  MOM  module  Is  the  third  module  run.  It  starts  up  from  the 
GTO  module  checkpoint,  then  begins  execution  as  shown. 
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GEMACS-MOM  TASK  EXECUTION  STARTED  ON  03/02/83  AT  10.03 

NUMBER  OF  PERIPHERAL  FILES  AVAILIBLE  17 

RUN  TIME  SET  TO  4.00  CPU  MINUTES 

FREQUENCY  SET  TO  500.  MEGAHERTZ 
HAVLENGTH  0.600  METERS 

FILL  IMPEDANCE  MATRIX  ZGARY 
USING  BASIS  FUNCTION  SINCOS 
ON  GEOMETRY  DATA  SET  THIELE 
SHADOWING  MATRIX  ZGSHD 

LOADS! IF  SPECIFIED) IN 
FREQUENCY (MEGAHERTZ )  500.00 

GROUND  COMO  (MHOS/M)  0. 

RELATIVE  PERMITIVITY  1.0000 

INTERACTIONS  SET:  PR  RR  PD  RO  PDR  PL  MM  EU  ES 

EXCITE  GEOMETRY  OATA  THIELE 
EXCITATION  VOLTGE 
EXCITATION  DATA  SRC 

INTERACTIONS  SET:  PR  RR  PO  RD  PDR  PL  MM  EU  ES 

REAL  COMP  1.00  IMAG  COMP  0. 

EXCITED  SEGS 
1 

DECOMPOSE  ZGARY  STORE  RESULT  IN  ZGARY  PIVOT-  N 

MAX  OIAG  -  15564.  MIN  DIAG  -  1580.9 

PIVOT  RATIO  -  9.85 

SOLVE  ZGARY*  I-  SRC 

GEOMETRY  OATA  SET  THIELE 

***  NO  LOAD  FOR  STRUCTURE  *** 

ANTENNA/LOAD  PARAMETERS  (SEE  USER  MANUAL) 

INPUT  INPUT 

SEGMENT  IMP(MAG)  IMP(PHZ)  PWR  INPUT  PWR  LOAD 
1  428.952  32.389  0.984E-03  0. 


SYMBOL  I 

LINEAGE-  SRC-  ZGARY-THIELE- 
COMPLX  DATA 


COLUMN-  1 

REAL  IMAGINARY  MAGNITUDE  PHASE(DEG)  REAL  IMAGINARY  MAGNITUDE  PHASE (DEG) 

1  0.1969E-Q2  -0.1249E-02  0.2331E-02  -32.39  2  0.1024E-02  -0.3394E-02  0.3545E-02  -73.22 

3  0.2492E-03  -0.4312E-02  0.4319E-02  -86.69  4  -0.2539E-03  -0.4235E-02  0.4242E-02  -93.43 

5  -0.4487E-03  -0.3310E-02  0.3341E-02  -97.72  6  -0.3037E-03  -0.1588E-02  0.1617E-02  -100.8 

GENERATING  FLOOIP  FROM:  INCIOENT  FIELD  MATRIX  FLDOIP 

GREENS  FUNCTION  MATRIX  FLDGFM 
SOLUTION  MATRIX  I 


LOADING  COLUMNS  1  TO  74  TO  GENERATE  COLUMN  1  OF  FIELO  MATRIX 


The  MOM  module  ends  execution  with  a  checkpoint,  from  which  the 
OUTPUT  module  begins  execution.  The  purpose  of  the  output  module  Is  to 
print  (and  plot)  the  field  patterns  generated  by  the  previous  modules. 

E-FIELO  MATRIX  FLDOIP 
SPHERICAL  COORDINATE  SYSTEM 

FAR  FIELO  FOR  FIELD  OATA-FLDOIP  -CURRENT/SOURCE  DATA-  I  -GEOMETRY  0ATA«THIELE 
NORMALIZATION  FACTOR  0.328  V/M 
CONSTANT  PHI-  0 

E(THETA)  El  PHI ) 


TH« 

MAGNITUDE 

PHASE (DEG) 

MAGNITUDE 

PHASE (DEG) 

NORMALIZED 

0. 

0. 

0. 

0. 

0. 

-100. 

5.0 

0.666E-01 

120. 

0. 

0. 

-13.8 

10.0 

0.103 

119. 

0. 

0. 

-10.0 

15.0 

0.937E-01 

114. 

0. 

0. 

-10.9 

20.0 

0.5Q2E-01 

87.2 

0. 

0. 

-16.3 

25.0 

0.621E-01 

2.86 

0. 

0. 

-14.4 

30.0 

0.132 

-14.1 

0. 

0. 

-7.89 

35.0 

0.192 

-12.1 

0. 

0. 

-4.63 

40.0 

0.234 

-4.30 

0. 

0. 

-2.93 

45.0 

0.263 

6.09 

0. 

0. 

-1.91 

50.0 

0.287 

16.7 

0. 

0. 

-1.16 

55.0 

0.307 

25.6 

0. 

0. 

-0.549 

60.0 

0.322 

31.7 

0. 

0. 

-0.141 

65.0 

0.328 

34.9 

0. 

0. 

0. 

70.0 

0.322 

35.5 

0. 

0. 

-0.152 

75.0 

0.306 

33.7 

0. 

0. 

-0.587 

80.0 

0.282 

29.7 

0. 

0. 

-1.29 

85.0 

0.253 

23.6 

0. 

0. 

-2.24 

90.0 

0.260 

0.202 

0. 

0. 

-2.00 

95.0 

0.240 

-10.2 

0. 

0. 

-2.71 

100.0 

0.221 

-21.4 

0. 

0. 

-3.40 

105.0 

0.203 

-32.9 

0. 

0. 

-4.14 

110.0 

0.181 

-44.6 

0. 

0. 

-5.13 

115.0 

0.153 

-57.0 

0. 

0. 

-6.63 

120.0 

0.116 

-72.7 

0. 

0. 

-9.03 

125.0 

0.762E-01 

-99.7 

0. 

0. 

-12.7 

130.0 

0.597E-01 

-154. 

0. 

0. 

-14.8 

135.0 

0.848E-01 

163. 

0. 

0. 

-11.7 

140.0 

0.112 

143. 

0. 

0. 

-9.33 

145.0 

0.115 

132. 

0. 

0. 

-9.12 

150.0 

0.848E-01 

121. 

0. 

0. 

-11.7 

155.0 

0.316E-01 

90.3 

0. 

0. 

-20.3 

160.0 

0.467E-01 

-27.7 

0. 

0. 

-16.9 

165.0 

0.928E-01 

-43.8 

0. 

0. 

-11.0 

170.0 

0.101 

-48.5 

0. 

0. 

-10.2 

175.0 

0.654E-01 

-50.5 

0. 

0. 

-14.0 

180.0 

0. 

0. 

0. 

0. 

-100. 

The  OUTPUT  module  concludes  with  a  checkpoint.  Should  the  user  wish 
further  results  from  this  problem  (another  pattern  cut  perhaps)  he  may 
begin  execution  from  the  OUTPUT  module  checkpoint  and  not  have  to  rework 
75  percent  of  the  problem. 


CHECKPOINT  NUMBER  4  STARTED  AT  TIME  0.079  ON  LOGICAL  UNIT 

COMMON  BLOCK  AOEBUG  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  AMPZIJ  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  ARGCOM  WRITTEN  OUT  TO  ICKFIL 

COMNM  BLOCK  CSV STM  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  DEFDAT  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  FLDCOM  WRITTEN  OUT  TO  ICKFIL 

COMNM  BLOCK  GEOOAT  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  IOFLES  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  JUNCOM  WRITTEN  OUT  TO  ICKFIL 

COMNM  BLOCK  PARTAB  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  SCNPAR  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  SEGMNT  WRITTEN  OUT  TO  ICKFIL 

COMNM  BLOCK  SYMSTR  WRITTEN  OUT  TO  ICKFIL 

COMMON  BLOCK  SVSFIL  WRITTEN  OUT  TO  ICKFIL 

COMNM  BLOCK  TEMPO!  WRITTEN  OUT  TO  ICKFIL 

COMNM  BLOCK  GTOOAT  WRITTEN  OUT  TO  ICKFIL 

COMNM  BLOCK  MODULE  WRITTEN  OUT  TO  ICKFIL 

COMNM  BLOCK  INTMAT  WRITTEN  OUT  TO  ICKFIL 

PERIPHERAL  FILE  8  SYMBOL  THIELE  NUMREC-  1 

PERIPHERAL  FILE  9  SYMBOL  ZGARY  NUMREC-  6 

PERIPHERAL  FILE  11  SYMBOL  SRC  NUMREC-  1 

PERIPHERAL  FILE  16  SYMBOL  I  NUMREC-  1 

PERIPHERAL  FILE  13  SYMBOL  FLDOIP  NUMREC-  1 

PERIPHERAL  FILE  10  SYMBOL  ZGSHD  NUMREC-  1 

PERIPHERAL  FILE  12  SYMBOL  FLOGFM  NUMREC-  74 

PERIPHERAL  FILE  14  SYMBOL  ZGUPR  NUMREC-  6 

PERIPHERAL  FILE  IS  SYMBOL  ZGLWR  NUMREC-  6 

CHECKPOINT  COMPLETE  AT  0.092 

ELAPSED  TIME-  0.012 

CURRENT  FILE  LENGTH-  25149 

STOP  AT  LINE  1 . 
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comwnications,  electromagnetic  guidance  and  control,  aut- 
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physics  and  electronic  reliability,  maintainability  and 
compatibility. 
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